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The driven transport of plastic systems in various disordered backgrounds is studied within mean 
field theory. Plasticity is modeled using non-convex interparticle potentials that allow for phase 
slips. This theory most naturally describes sliding charge density waves; other applications include 
flow of colloidal particles or driven magnetic flux vortices in disordered backgrounds. The phase 
diagrams exhibit generic phases and phase boundaries, though the shapes of the phase boundaries 
depend on the shape of the disorder potential. The phases are distinguished by their velocity and 
coherence: the moving phase generically has finite coherence, while pinned states can be coherent or 
incoherent. The coherent and incoherent static phases can coexist in parameter space, in contrast 
with previous results for exactly sinusoidal pinning potentials. Transitions between the moving and 
static states can also be hysteretic. The depinning transition from the static to sliding states can 
be determined analytically, while the repinning transition from the moving to the pinned phases is 
computed by direct simulation. 

PACS numbers: 83.60.Bc,62.20.Fe 



I. INTRODUCTION 

The collective dynamics of extended systems driven 
through quenched disorder is a rich and challenging prob- 
lem, with many experimental realizations. Such systems 
include vortices in type II superconductors, charge den- 
sity waves in anisotropic conductors, domain walls in ran- 
dom ferromagnets, and planar cracks in heterogeneous 
materialsii Much of the theoretical work to date has fo- 
cused on modeling these systems as extended elastic me- 
dia. In these models the restoring forces are monoton- 
ically increasing functions of the relative displacements, 
and the system is not allowed to tear. At zero temper- 
ature, overdamped elastic media subject to an applied 
force F and quenched disorder exhibit a nonequilibrium 
phase transition from a pinned state to a sliding state 
at a critical value, Ft, of the driving forcei^ The depin- 
ning transition, first fully studied for collective models 
with disorder in the context of charge density waves, dis- 
plays the universal critical behavior of continuous equi- 
librium phase transitions, with the mean velocity v of the 
medium playing the role of the order parameteri^ For 
monotonic interactions, it has been shown that the sys- 
tem's velocity is a unique function of the driving forceA 
The sliding state is therefore unique and there is no hys- 
teresis or history dependence. The depinning transition 
of driven elastic media has been studied extensively, both 
by functional renormalization group methods^*^*^ and 
large scale numerical simulationsiS*2*i2iii4i2ii^ Universal- 
ity classes have been identified, which are distinguished, 
for example, by the range of the interactions or by the pe- 
riodicity (or nonperiodicity) of the pinning force. More 
recent work, while still focusing on elastic media, has 
shown that the dynamics is quite rich well into the uni- 
formly sliding state. "•^^■^'^•^'■18.20 

The elastic medium model is often inadequate to de- 
scribe many real systems which exhibit plasticity (due, 
for instance, to topological defects in the medium) or in- 



ertial effects that violate the assumption of overdamped 
equations of motion. The dynamics of plastic systems can 
be both spatially and temporally inhomogeneous, with 
coexisting pinned and moving regions d^ The depinning 
transition may become discontinuous (first order), possi- 
bly with macroscopic hysteresis and "switching" between 
pinned and sliding statesiSiS^iS^iSi The theoretical un- 
derstanding of the dynamics of such "plastic" systems is 
much less developed than that of driven elastic media. A 
number of mean-field models of driven extended systems 
with locally underdamped relaxation or phase slips have 
been proposed in the literature)^i^^i^^i^^i^^i^^i^°i^^i^^i'^^ 
but many open questions remain. 

Much of the original theoretical work on driven dis- 
ordered systems was motived by charge density wave 
(CDW) transport in anisotropic conductors, which dis- 
play a nonlinear current-voltage characteristic with a 
threshold voltage for collective charge transport. ■^^■'^^ It 
has been known for some time that the elastic depin- 
ning transition may not be physically relevant to real 
CDW materialsi^^*^Si2i Coppersmith argued that in elas- 
tic models with weak disorder, unbounded strains can 
build up at the boundaries of an atypically low pinning 
region, resulting in large gradients of displacement that 
lead to the breakdown of the elastic model."^^ Topolog- 
ical defects or phase slips will occur at the boundaries 
of such a region, yielding a spatially nonuniform time- 
averaged velocity. Theoretical and numerical studies of 
models that incorporate both phase and amplitude fluc- 
tuations of the CDW order parameter have indicated 
that phase slips from large amplitude fluctuations can 
destroy the critical behavior iSSiSSi^ The depinning may 
become discontinuous and hysteretic, or rounded, in the 
inflnite system limit. Experiments show that varying 
the temperature of the CDW material can lead to a 
transition from continuous depinning to hysteretic de- 
pinning with sharp "switching" between pinned and slid- 
ing statesiSSiiSi^ Furthermore, the observed correlation 
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between the amplitude of broadband noise and macro- 
scopic velocity inhomogeneities also suggest the presence 
of phase slipsi^ It should be mentioned, however, that 
in many samples a substantial amount of phase slips oc- 
curs at the contacts,'*'^ while less clear evidence exists 
for substantial phase slip effects in the bulk. In gen- 
eral, CDW experiments display considerable sample-to- 
sample variability,^^ making the comparison between the- 
oretical models and experiments quite challenging. 

Related slip effects or plastic behavior have been 

proposed to explain the complex dynamics of many 

other dissipative systems, including vortex arrays in 

type-II superconductors. Simulations (mainly in two 
dimensions) I7i44i45,46,47,48,49^ iniagingi8-5"-5i^52^, and 

transport and noise experiments^^'^S'^ have shown that 
driven flux lattices often do not respond as elastic me- 
dia. Instead, the driven lattice tears as small-scale topo- 
logical defect structures are generated and healed by 
the interplay of drive, disorder and interactions. The 
tearing results in a "plastic" response, with highly de- 
fective liquid-like regions flowing around the boundaries 
of pinned solid-like regions."*^ This kind of response is 
most prominent in the region near vortex lattice melt- 
ing, where the so-called peak effect occurs, i.e., the crit- 
ical current shows a sudden increase with temperature 
or applied field. Reproducible noise or "fingerprint phe- 
nomena" have been observed in the current-dependent 
differential resistance and attributed to the sequential de- 
pinning of various chunks of the vortex lattice,^"* Images 
of driven vortex arrays in irradiated thin films of Niobium 
obtained by Lorentz microscopy have shown clearly that 
vortex rivers flowing past each other at the boundaries 
of pinned regions of the lattice. Scanning tunneling mi- 
croscopy, which can resolve individual vortices at high 
density, has also revealed a clear evolution of the vor- 
tex dynamics with disorder strength,-^^ In samples with 
weak disorder the vortex array was observed to creep co- 
herently along one of the principal crystal axes near the 
onset of motion. In samples with strong disorder, the de- 
pinning is plastic and the path of individual vortices can 
be followed as they meander through the pinned crystal. 
Finally, as in the case of CDWs, a correlation between 
plasticity and broadband noise has been observed in sev- 
eral samples. Recently it has been argued that some 
of the observed behavior may be due to edge contam- 
ination effects that are responsible for the coexistence 
of a metastable disordered phase and a stable ordered 
p]^g^gg|57j58j59 jg Qiga,r that more work is needed to un- 
derstand the rich dynamics of these driven systems. 

In this paper we study the driven dynamics of a dis- 
ordered medium with phase slips, in order to better ad- 
dress questions about these and related physical systems. 
We restrict ourselves to systems which are periodic along 
the direction of motion, such as CDWs, vortex lattices or 
2D colloids, and consider only the dynamics of a scalar 
displacement field. For concreteness, the model is de- 
scribed in the context of driven CDWs, but it also ap- 
plies to other driven systems with pinning periodic in 



the displacement coordinate. Assuming overdamped dy- 
namics and discretizing spatial coordinates, the dynam- 
ics of the phase 9i of each CDW domain is controlled 
by the competing effects of the external driving force, 
the periodic pinning from quenched disorder, and the 
interaction among neighboring domains. Following the 
literaturer^iSfliSiiSS phase slips are introduced by model- 
ing the interactions as a nonlinear sine coupling in the 
phase difference of neighboring domains. The mean field 
limit for this type of model has been studied by Stro- 
gatz, Wester velt, Marcus and MiroUo^^ for the case of 
the smooth sinusoidal pinning force and was shown to 
exhibit a first order depinning transition, with hysteresis 
and switching. In this paper we use a combination of ana- 
lytical methods and numerical simulations to obtain the 
nonequilibrium mean field phase diagram of the phase 
slip model for a variety of pinning forces (see Fig. 
Note that most of the pinning forces we consider are dis- 
continuous. This form of the force mimics the cusped 
potentials that are the starting points for mean field the- 
ories that best reproduce the finite-dimensional results. 
The discontinuous pinning forces also reflect the abrupt 
changes in the effective force (sum of elastic and pin- 
ning forces) that occur when a neighboring region of the 
medium suddenly moves forward. We find that discontin- 
uous forces, and even continuous nonsinusoidal pinning 
forces, yield a rich nonequilibrium phase diagram, with 
novel stable static phases that are not present for exactly 
sinusoidal pinning forces. 

In mean field theory, the nonequilibrium state of the 
system can be described in terms of two order parame- 
ters. As the pinning potential for each domain i is pe- 
riodic in 9i, having minima at [3i + 27rn, for integer n, 
and taking the interactions to be periodic in the differ- 
ence 9i — 6j between neighboring phases with the same 
period, a natural order parameter is the coherence of the 
phases. This coherence is measured by the amplitude r 
of a complex order parameter defined via 

1 ^ 

with ip a mean phase. In the absence of interactions 
among the phases or external drive, the 9i's are locked to 
the random phases, 9i = Pi, and the state is incoherent, 
with r = 0. In the opposite limit of very strong interac- 
tions we expect perfect coherence of the static state, with 
all phases becoming equal and r — > 1 as the interactions 
become strong (or the pinning becomes weak.) Another 
order parameter is the average velocity of the system, 
given by 

1 ^ 

^=]^E^^(^)- (2) 

The mean velocity is the order parameter for the transi- 
tion between static and moving phases. 
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FIG. 1: Sketches of the pinning potentials and forces studied 
in this paper. The pinning forces are periodic with period 
2ti and the pinning potential for a degree of freedom 9i has 
minima at Pi + 2n7r, for integer n. The cases are organized 
primarily by the sign of c, with the pinning force Y(x) — 
—ax — cx^ + 0{x^) for small x = 6i — f3i. The coefficient 
of the harmonic part of the force satisfies a > 0. The cases 
(a), (b) and (c) are for "soft" pinning forces (c < 0); they 
differ near the potential maxima, corresponding to monotonic, 
nonmonotonic, and continuous forces, respectively. Case (d) 
is a "hard" potential (c > 0). The "scalloped" potential, case 
(e), is precisely quadratic (c = 0) in the interval — vr < a; < tt. 
The form of the potential especially affects the stability of the 
coherently pinned phase and whether "reentrant" pinning is 
possible upon increasing or decreasing the force. 



The central results of this paper are the nonequilib- 
rium phase diagrams describing the static and moving 
phases, for the various pinning forces shown in Fig. ^ 
The parameters for the phase diagrams are the driving 
force F and the strength /i of the interaction between the 
domains. (For a phase diagram in the drive force vs. pin- 
ning strength plane, see Sec lVIIl l Although the precise 
shape of the phase boundaries depends on the detailed 
form of the pinning potential, the types of phases and 
the schematic topology of the phase diagram are gen- 
eral. This topology and set of phases is exemplified in 
the phase diagram for the discontinuous soft cubic pin- 
ning force (see Fig.^b)) shown in Fig.|21 We find three 
distinct zero-temperature nonequilibrium phases: 

• an incoherent static phase (IS) at low drives and 
small coupling strengths, with v = and r — 0; 

• a coherent static phase (CS) at low drives and large 
coupling strengths, with v = and r > 0; 

• a coherent moving phase (CM) at large drives, with 
V > and r > 0. 

We have investigated the possibility of an incoherent 
moving (IM) phase. For continuous pinning forces, there 
is no IM phase. For discontinuous pinning forces, we 
speculate that the IM phase is unstable generically. (See 




FIG. 2: Phase diagram in the coupling-drive (i-i-F) plane for 
a discontinuous soft cubic pinning force of the type shown 
in Fig. The equation of motion is Eq. The corre- 

sponding Y{x) is given by Eq. 13411 with a = 15/(87r) and 
c = — 4a^/27. The strength of the pinning is h — 1 for all de- 
grees of freedom. The diagonally lined region indicates the IS 
phase, while the cross-hatched region indicates the CS phase. 
The light gray shaded region denotes the region of coexistence 
of the CM and IS phases, while the medium gray shaded re- 
gion denotes the region of coexistence of the IS and CS phases. 
The lines and F^ are the forces at which the system de- 
pins upon increasing the drive from the incoherent and co- 
herent static states respectively. The line F| is the force at 
which a coherently moving system stops upon lowering the 
drive. The point {fie,Fe) indicates where the static-moving 
transition goes from hysteretic to non-hysteretic. The curves 
fiu{F) and iJ.d{F) are the values of the coupling at which the 
static system makes the transition to and from finite coher- 
ence states, respectively. The inset displays the hysteresis in 
the coherence r as the coupling strength jj, is varied at _F = 0. 
The transitions between the IS and CS phases are first order 
in r. 



Sec. V where the stability of a possible IM phase is dis- 
cussed.) 

An important new feature of the phase diagram is the 
occurrence of a coherent static phase at finite F. In con- 
trast, for the sinusoidal pinning force studied previously 
by Strogatz and collaborators^^ the static state is always 
incoherent (IS) for all finite values of the driving force 
and the CS phase is only present at F = 0. 

The location of the transitions between these phases 
depends on the system's history. Changing the coupling 
fi at fixed drive F can give a hysteretic transition be- 
tween incoherent and coherent static phases, as shown in 
the inset of Fig. |21 for F = 0. Fig. |3| shows the behav- 
ior of both the mean velocity and the coherence as F is 
first increased and then decreased across the boundaries 
between static and moving phases of Fig. |21 while keep- 
ing fi fixed. The most important features of the phase 
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diagrams are: 



• The transition between the IS and CS phases is gen- 
erally discontinuous. The region of coexistence of 
coherent and incoherent static states is bounded 
by curves ^id{F) and ^iuiF) (or equivalently Fd{^) 
and F^ip)-) When the couphng strength ^ is in- 
creased at fixed F within the static region, the sys- 
tem jumps from an incoherent to a coherent state 
at the critical value fiu{F), with a discontinuous 
change in r (see inset of Fig.|21l. When ^ is ramped 
back down, the coherent static state remains stable 
down to the lower value fid{F). The boundaries 
IJ,d{F) and /i«(F) coincide for the piecewise linear 
pinning force. In this case the transition is still 
discontinuous, but not hysteretic. An exception to 
this general behavior is found for the hard pinning 
potential at very small values of F, where the tran- 
sition between coherent and incoherent static states 
is continuous. 

• The depinning to the moving phase is discontinu- 
ous and hysteretic when the system depins from the 
IS phase (except when ^ — 0). When F is increased 
adiabatically from zero at fixed fi for a system pre- 
pared in the IS phase, both the velocity and the 
coherence jump discontinuously from zero to a fi- 
nite value at FUji). For an example, see the top 
frames of Fig. y] When the force is ramped back 
down from the sliding state the system gets stuck 
again at the lower value F^ (/i) . 

• The depinning to the moving phase is generally 
continuous when the system depins from the CS 
phase. In this case both the velocity and the co- 
herence change continuously at the transition, al- 
though they may be non-analytic functions of the 
control parameters. An example of this behavior is 
displayed in the bottom frames of Fig. |21 An ex- 
ception is found for piecewise linear pinning forces 
(case (e) of Fig. for fi^ fiu- 

• For continuous pinning forces, the depinning 
threshold F^(fj,) vanishes for fi above a critical /iy. 
In contrast, discontinuous pinning forces exhibit a 
finite depinning threshold for all finite values of /i 
with F^(fj,) decreasing with increasing fi. 

Analytical expressions have been obtained for the criti- 
cal lines F^{^) and F^{fi), which give the depinning force 
values for the coherent and the incoherent static phases, 
respectively, as well as for the phase boundaries fJ-d{F) 
and fiu{F), which separate the coherent and incoherent 
static phases. Numerical simulations of finite mean-field 
systems have also been used to obtain these boundaries, 
confirming the analytic stability criteria. The repinning 
curves (i^|(yLt)), where moving solutions stop upon lower- 
ing the drive F, have been determined numerically. 

Part of the motivation for our work comes from the 
well-known result that the mean field critical exponents 
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FIG. 3; Typical numerical results, found by integrating nu- 
merically the equations of motion (Eq. l|3J), for the behavior 
of the mean velocity v and the coherence r as the driving 
force is slowly varied. For each pair of plots, the coupling fi 
is held constant, while the drive force F is raised from F = 
to F = 1.2 and then decreased. The pinning potential is the 
same as for Fig. |21 The top frames (/x = 0.5) show the hys- 
teretic behavior between the IS and the CM phases, where the 
coherence and velocity jump between zero and non-zero values 
at the same locations. The next two sets of frames — 1.14) 
are obtained by preparing the system in the IS-CS coexistence 
region, starting from either an initial incoherent (/) or coher- 
ent (C) state. When the system is prepared in an incoherent 
state, the velocity and coherence jump at at the same value 
of F (« 0.42) as F is raised, but change continuously as F 
is decreased, albeit with a change in the slope dr/dF at the 
repinning force F ~ 0.32, where v goes to zero. When the sys- 
tem is prepared in a coherent state, there is no hysteresis and 
V and r are continuous, though r again shows a singularity at 
depinning. The bottom frames = 1.5) display the behav- 
ior at the continuous depinning transition from the CS phase. 
The results are similar to those for = 1.14, when starting 
from the coherent state (C). In general, depinning from the 
coherent state is continuous and non-hysteretic, while depin- 
ning from the incoherent state is discontinuous and hysteretic. 
Numerical evidence for the hysteresis does not change over the 
size ranges studied, strongly suggesting that these simulations 
accurately represent the infinite-volume limit. 
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for the depinning transition in purely elastic models de- 
pend on the details of the pinning force. For instance, the 
exponent (3 controlling the vanishing of the mean velocity 
V with driving force at threshold, w ~ (F — Ft)^ , has a 
mean field value (3 — 3/2 for generic smooth continuous 
pinning forces and /3 = 1 for a discontinuous piecewise 
linear pinning force (Fig. 1(e)). Using a functional RG 
(FRG) expansion in 4— e dimensions, Narayan and Fisher 
showed'^ that the discontinuous force captures a crucial 
intrinsic discontinuity of the large scale, low-frequency 
dynamics. The FRG calculations give /3 = l—e/6+0{e'^), 
in good agreement with numerical studies in two and 
three dimensions. The mean field elastic medium also has 
zero depinning field. Ft — 0, for small pinning strengths 
h, in contrast with finite dimensional simulations and 
predictions for a finite depinning field in any dimension 
based on Imry-Ma/Larkin-Ovchinnikov and rare region 
argumentsj2. The RG calculation and the numerics show 
that a discontinuous pinning force must be used in the 
mean field theory to incorporate the inherent jerkiness of 
the motion of finite-dimensional systems at slow veloci- 
ties. Although there is no reason to believe a priori that 
the same will hold for models with phase slips, it is clearly 
important to understand how the properties of the pin- 
ning potential affect the nonequilibrium phase diagram 
of the model. Furthermore, for large coupling strength fi 
and bounded pinning force the phase slip model reduces 
to the elastic model, where the nature of the pinning force 
strongly affects the mean field theory. 

For further applications and connections, we note that 
models of driven disordered systems with nonmonotonic 
interactions are also relevant for arrays of nonlinearly 
coupled oscillators. An example is the Kuramoto model 
used to describe the onset of synchronization in many bi- 
ological and chemical systems.^- The model consists of a 
large number of oscillators with random natural frequen- 
cies and a sinusoidal coupling in their local phase difi'er- 
ences. Although there is no external drive, this model 
can exhibit a transition to a synchronized phase as the 
strength of the coupling is increased. In this phase, all 
the degrees of freedom oscillate at a common frequency. 
In the Kuramoto model the natural frequency acts as a 
random driving force that varies for each oscillator, but 
there is no random pinning. The model considered here, 
in contrast, consists of coupled phases, or oscillators, in 
a random pinning environment at fixed (constant) drive. 
The onset of coherence (either in a moving or in a static 
state) corresponds to the onset of the synchronization in 
the Kuramoto model. 

We conclude this introduction by briefly summarizing 
the remainder of the paper. In Sec. ^1 we describe the 
model of driven CDWs with phase slips and introduce 
the mean field limit. In Sec. IIIII we obtain the static so- 
lutions of the mean field model at F = for the selection 
of pinning forces shown in Fig. We show that the ex- 
istence of a transition between incoherent and coherent 
static states can be inferred perturbatively. A full non- 
perturbative treatment is then applied to understand the 



nature of the transition. In Sec. IIVI we consider static 
states at finite drive. Again, the region of stability of the 
incoherent static phase can be established by perturba- 
tion theory, but the nonperturbative treatment described 
in Sec. is needed to map out all the static states and 
their boundaries of stability to the moving state. The re- 
sulting phase diagrams for the various classes of pinning 
forces are discussed in Sec. the analytic calculations 
supporting these phase diagrams are presented in Ap- 
pendices El ^nd^ As the analytic treatment we present 
here is restricted to finding boundaries starting from the 
static phases, the lower boundary Fx(/^) of the hysteretic 
region where static and moving state coexist has been 
obtained numerically. Sec. IVII addresses the effect of a 
broad distribution of pinning strengths. We conclude in 
Sec. IVIII with a discussion of the results and avenues for 
further studies. 



II. THE MODEL 

Though the results of our analysis are more general, 
we motivate the model with a detailed discussion of the 
physics of CDWs. The general ideas of phase slip also ap- 
ply to other systems, most directly to coupled layers of 
vortices, where the vortices are confined to the planar lay- 
ers, or to colloidal particles in a disordered background. 

A CDW is a coupled periodic modulation of the elec- 
tronic density and lattice ion positions that exists in cer- 
tain quasi-one dimensional conductors, due to an insta- 
bility of the Fermi surface. The undistorted CDW state 
is a periodic condensate of electrons, characterized by 
a complex order parameter, with an amplitude pi and 
a phase 9. The electron density can be expanded as 
Pe{x) = po + piCOs[QcX + 0{x)], with Qc = 2kp, kp 
being the Fermi wave vector. The phase 9{x) describes 
the position of the CDW with respect to the lattice ions 
and is a constant for an undistorted CDW. When Qc is 
incommensurate with the lattice, the CDW can "slide" 
and CDW transport can be modeled using uniform trans- 
lations and small gradients of 9{x), to a first approxi- 
mation. An applied electric field exceeding a threshold 
field causes the CDW to slide relative to the lattice at 
a rate dt9, giving rise to a CDW current. Amplitude 
fluctuations (changes in pi) are often neglected because 
they cost a finite energy, while a vanishingly small en- 
ergy is required to generate long-wavelength phase ex- 
citations, in an ideal crystal. This has led to the well 
known phase-only model of CDW dynamics introduced 
by Fukuyama, Lee and Rice (FLR) that incorporates long 
wavelength elastic distortions of the phase. ®^ Strong dis- 
order or regions of unusually low pinning can lead to large 
strains, however, so that the amplitude can no longer be 
regarded as constant. Large local strains can be relieved 
by a transient collapse of the CDW amplitude. One ap- 
proach to describe such a strongly strained system is a 
"soft spin" model that considers the coupled dynamics 
of both phase and amplitude fluctuations. This has been 
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attempted by some authors r2i^2i2S but generally leads to 
models that have to be treated numerically. An alterna- 
tive, more tractable approach, is to continue to treat the 
amplitude as constant, while modifying the interaction 
between phases. This modification should incorporate 
the crucial feature that the phase becomes undefined at 
the location where the amplitude collapses. At a strong 
pinning center, phase distortions can be large and lead to 
the accumulation of a large polarization that suppresses 
the CDW amplitude, driving it toward collapse. When 
the distortion is released through an amplitude collapse, 
the phase abruptly advances of order « 27r, while the 
amplitude quickly regenerates*^ This process is known 
as phase slippage in superconductors and superfluids, al- 
though it is modified in CDWs because of the physical 
coupling to the phase. On time scales large compared to 
those of the microscopic dynamics, it can be described 
approximately as a "phase slip" : an instantaneous 27r 
(modulo 27r) hop of the CDW phase. Following the liter- 
ature, phase slips will be modeled here as phase couplings 
periodic in the phase difference between neighboring do- 
mains. This leads to a simple model that can be analyzed 
in some detail. 

When modeling CDWs, especially numerically, dis- 
placements and amplitudes are coarse-grained to a length 
scale of order the Imry-Ma-Larkin-Ovchinikov length. At 
and below this scale, the CDW behaves roughly as a rigid 
object, referred to as a correlated domain. This domain 
is taken to move uniformly and is acted upon by driving 
forces and interactions with neighboring domains and the 
pinning potential. The continuum space description is 
replaced with a discrete set of degrees of freedom. The 
coarse-grained equation of motion for the phase 9i of a 
CDW domain i is given by 

e, = F + ^lY^ sin(0j - Oi) + h,Y{e, -(3i) , (3) 

where the dot denotes the time derivative (we have cho- 
sen to scale time so that the damping constant is unity) 
and F is the driving force. The second term on the right 
hand side of Eq. Q represents the force due to the cou- 
pling to other domains, where (j) ranges over sites j that 
are nearest neighbor to i and /i is the coupling strength. 
The third term is the pinning force which tends to pin 
the phase of each domain to a random value /3i uniformly 
distributed in [— tt, tt]. The function Y{x) is periodic with 
period 27r and represents the pinning forces. We choose 
y(0) = to fix the location of the minimum of the pin- 
ning potential and set F"(0) = to maintain reflection 
symmetry in the absence of an external drive. As the 
potential is minimized at Oi — Pi, Y'{Q) < 0. The ran- 
dom pinning strengths hi are independently chosen from 
a probability distribution p{h). 

The key difference between our model equation of mo- 
tion and the well known FLR elastic model of driven 
CDWs is in the form of the coupling between domains. 
Instead of assuming a linear elastic force ^ Tli<j> 
between neighboring domains, we have assumed a non- 



linear, sinusoidal coupling that allows for phase slip pro- 
cesses. For large phase distortions (exceeding tt) the 
restoring force in Eq. Q becomes negative and the 
phases slip by an amount 27r relative to one another in 
order to relax the strain. 

The starting point for many finite-dimensional theo- 
ries is the mean field picture where every local phase (or 
domain) is equally coupled to every other. In this limit, 
the equation of motion (jSJ becomes 

eiP, h)^F -u sm{e - V') + hY{e - P) , (4) 

where 

u = nr , (5) 

measures the effective strength of coupling between the 
domains and the mean field, with r and ip defined in 
Eq. £3). This coupling will only be non-zero if there is 
some coherence between the phases of different domains, 
i.e., if r ^ 0. For simplicity, we have dropped the sub- 
scripts, labeling each phase by the values of /? and h, 
which are now both continuous variables. The /? are dis- 
tributed uniformly in [—tt, tt] and the h have the distri- 
bution p{h). 

The self-consistency condition for the mean field theory 
is given by 

re'^ f dp j d/ip(/i)e'^('3,/.) (g) 

In this paper we will for the most part consider a narrow 
distribution of pinning strengths, i.e., p(h) = 5{h — 1). 
The effects of a broad distribution p{h) will be addressed 
in Sec. ED 

When the phases are not coupled [p — 0), the equa- 
tion of motion reduces to that of a single particle, which 
depins at the single particle threshold force, Fgp, given 
by the maximum pinning force. Note that when the co- 
herence r is zero, then u = 0, and the system may also 
depin at _Fsp for a finite value of p, as long as r remains 
zero. 



III. STATIC STATES FOR ZERO DRIVE 

We first consider static solutions {9 = 0) to Eq. Q 
for the case of zero drive {F = 0). These solutions are 
the first step in determining the phase diagram and their 
derivation introduces most of the techniques and con- 
cepts used for non-zero drive. When _F = 0, the coher- 
ence r is determined by competition between two effects: 
the disordering effect of the random impurities and the 
ordering tendency arising from the coupling of each de- 
gree of freedom to the mean field. The outcome of this 
competition gives the /^-dependence of r. At zero drive, 
the system can exist in one of two possible phases: the 
disordered (r = 0) IS phase and the ordered (r > 0) 
CS phase. These phases can coexist. In this section we 
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examine the nature of the transition between these two 
phases obtained by varying ^ at F = 0. We find that 
the nature of the transition depends on the shape of the 
pinning force, Y(x). 

For static solutions at zero drive, the equation of mo- 
tion 0) reduces to the condition that the pinning force 
on each degree of freedom be balanced by the force due 
to deformations from coupling to the mean field, 

0= -u sm{0 -TP) + hY{9~(3) , (7) 

where the reader is reminded that the effective coupling 
u results from the coupling strength fi and coherence r, 
u — iir. For any value of /i this equation has the trivial 
solution = P, r = u = 0, where all phases rest at the 
minima of their pinning potentials and the coherence and 
effective coupling are both zero. It turns out, however, 
that such a static incoherent solution becomes unstable 
above a characteristic value of the coupling strength /i. 

In order to study the competition between the impurity 
disordering and mean- field ordering effects, it is useful to 
rewrite the equation in terms of the deviation S of each 
phase from its value in the disorder dominated incoherent 
state, 5 = 6 — (3. A direct and important symmetry of the 
solution of Eq. Q is global phase invariance, which holds 
due to the uniform choice of /3. In the static state, this 
statistical rotational invariance means that we can simply 
fix -0 to be zero. Given a solution with = 0, all related 
solutions with ^ can then be obtained by letting 
6 ^ 9 + ip and (3 ^ (3 — i>. With this transformations, 
and specializing to the case of fixed pinning strength, 
h ~ 1, the force balance equation becomes 

= -usm{6 + 13) +Y{S) . (8) 

To solve this force balance equation, we need to deter- 
mine u self-consistently. The self-consistency condition 
Eq. © can be rewritten, by separating out its real and 
imaginary parts, as 

r=^ f dPcos{S + p)=f{u) , (9) 

where we have implicitly used Eq. JHJ to solve for S as 
a (possibly multi-valued) function of (3 and u to define a 
function f{u) as the above average over /?, and 

0=1 df3sm{S + f3) . (10) 

Next, we will use a straightforward linear stability anal- 
ysis to show that the IS (r = 0) phase becomes unstable 
to the CS (r > 0) phase above a critical value /i„ of the 
coupling strength. A perturbative calculation of r(fi) al- 
lows us to establish that this transition from the IS to 
the CS phase is continuous or hysteretic, depending on 
the shape of the pinning potential near its minimum. We 
will then obtain the full solution 7'(/i, _F = 0) for a variety 
of pinning forces. 



A. Stability of the Incoherent Static Phase 

To investigate the linear stability of the IS phase, we 
calculate the time evolution of a configuration near the 
static solution S{/3) = 0. A convenient perturbed config- 
uration is S{f3,t = 0) = -e(0)sin/3 with e(0) < 1. This 
perturbation gives non-zero coherence while maintaining 
-0 = and refiects the most rapidly growing eigenvector 
in the stability analysis, with S{(3, t) — — e(t) sin/3 to low- 
est order in e. By Eq. (O, the coherence of the perturbed 
state is 

1 r 

r ~ — / d/3cos(/3 — esin/?) , 

^e/2 + 0{e'). (11) 
The equations of motion Eq. Q) then give 

d = - ^ir sm{(3 + 6) + hY (6) 

= -fiie/2)sm(3 + hY'{0)S + 0{e'^) 

= {^ + hY'{0))5 + O{e'). (12) 

As r and S are both proportional to e (to lowest order), it 
immediately follows that f -|- /iF'(O)) r. The criti- 

cal value of 1-1 for linear stability is therefore 

flu = -2hY'{0) . (13) 

For coupling strength /j, > the perturbed coherence 
grows and the IS phase is linearly unstable to a CS phase. 
At larger /i, the interactions that drive the 9 towards a 
coherent state are larger in magnitude than the restoring 
force for the individual 9. Note that fi^ depends only on 
the strength of the pinning force at the minimum of the 
pinning potential. 



B. Perturbation Theory 

The onset of coherence for /i just above /i„ can be 
studied pcrturbatively by assuming that both the phase 
S and the coherence r are small in this region. Near 
(5 = the pinning force can quite generally be written as 
a power series in S, 

Y{S) = -aS -bS^ -cS^ + ... , (14) 

with a — ~Y'{0) > 0. For small r, and hence u, one can 
expand 6{(3,u) in powers of u, 

S{l3,u) ^ u6i{P) + u^52if3) + u'^SsiP) + ... . (15) 

Substituting these terms into the force balance equation 
Eq. ijHJl, and equating terms of the same order in u, we 
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obtain 



S3W 



sin/3 
a 

sin f3 cos(/3) b sin^ /3 

-i^ ] sin3/3 
sin/3cos^/3 36 sin^ /3 cos /? 



262 



(16a) 
(16b) 

(16c) 



Substituting the expanded 5(/3, m) into Eq. and eval- 
uating the integrals to each order in u we find 



f{u) = uri + u'^r2 + u^r^ + . . 



with 



1 

Ya ' 
r2 = 0, 

3 ( ac~2\? 



(17) 

(18a) 
(18b) 

(18c) 

Finally, the coherence r is given by the solution of 

r = /(r/i) = (r/i)ri + (r^) V3 + . . . . (19) 

For simplicity of discussion we specialize to pinning po- 
tentials with reflection symmetry and choose 6 = (al- 
though the non-zero 6 result will prove useful in the 
analogous finite F perturbation theory). Then = 
— 3c/(8a'*) and the nonvanishing solution for the coher- 
ence can be written as 




[Ai^] c<0, 

V {/2 (20) 
0, 



where ^.u = '2-a. 

The behavior of r(/i) for /i « /in and the nature of the 
transition between the IS and CS phases are controlled 
by the sign of the coefficient c of the cubic term oiY{6). 
The three types of behavior that can occur are shown in 
Fig. 0] For c > 0, corresponding to a "hard" pinning 
potential that grows more steeply than a parabola near 
its minimum, the coherence r grows monotonically with 
increasing /i, with r ~ (/x — ^uY^'^- This indicates a 
continuous transition at /i = between the IS and CS 
phases. On the other hand, when c < 0, corresponding to 
a "soft" pinning potential, the coherence starts out with 
a negative slope at /i„ and grows with decreasing /i. We 
expect this solution to be unstable, indicating that the 
transition from the IS phase to the CS phase occurs with 
a discontinuous jump in r from r = for ^ < to a 
non-zero value of r for fi > fi^ on a stable upper branch 
not accessible in perturbation theory. In fact we show 
below that when fi is decreased back down through /i„ 
from the CS phase r will remain non-zero down to a lower 



c = 




FIG. 4: The behavior of the coherence r for couplings /x « 
i.e., near the instability point of the incoherent static phase 
(IS) at _F = 0. The three curves show r with pinning force 
Y{5 ^ 1) = —aS — cS^ for c positive (hard pinning potential), 
negative (soft pinning potential) and zero (piecewise linear 
pinning force.) 



value fid < indicating a hysteretic transition between 
the IS and CS phases. In the marginal case of piecewise 
linear pinning forces with c = 0, i.e., Y{5) = —aS near 
S — 0, there is a discontinuous jump r{fi) at fi — fiu- In 
this case the perturbation theory breaks down and the 
solution must be obtained by the method described in 
Sec. nil CI This calculation will show that no hysteresis 
occurs in the case of strictly linear pinning force. We 
stress that the transition from the IS to the CS state at 
= is controlled entirely by the shape of the pinning 
potential near its minimum. Specifically, the behavior 
is unaffected by the existence of discontinuities in the 
pinning force at the edges of each pinning well. 



Beyond Perturbation Theory: The General 
Static r{fi, F = 0) Solution 



In this section we outline a non-perturbative method 
for calculating the integral f{u) used in the self- 
consistency equation, Eq. @. This allows for the de- 
termination of the coherence r for all values of fi. In 
addition to confirming the perturbative results obtained 
above, this method allows the precise study of the discon- 
tinuous and hysteretic transitions between the IS and CS 
phase, which cannot be done within perturbation theory. 

To obtain f{u) by direct integration over /3 in Eq. 
one would need to solve the transcendental equation, 
Eq. ijHJl, for S{P,u). Such a solution cannot in general 
be obtained analytically. Hence we take an alternative 
approach in which we solve Eq. ijSJ) for /3((5, u) and inte- 
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grate over S, rather than /3, i.e., 

r^^Jds(^^yos{5 + P(S,u)). (21) 

The change of variable in Eq. (|21|l provides an important 
simphfication that allows us to calculate analytically the 
coherence of the undriven static state for a general pin- 
ning potential. This simplification does rely on under- 
standing the subtleties of how S depends on (3, as S can 
be multivalued function of /3. The history of the sam- 
ple can determine which branch(es) are included in the 
configuration. 

For a given w, there is an infinite set of solutions to 
Eq. (jSJ. We index each with an integer n: 

/3n{6,u) = -5 + rnr +{-!)'"■ fim-\Y{d)/u) , (22) 

where we choose the [— 7r/2,7r/2] branch for sin^^(a;). 
The range for S is constrained to — (5,„ax(i*) ^ S < 

<5max(w), with (5max(M) = Y~'^{u). 

The calculation of the average in Eq. H21(l is easily car- 
ried out when u < a, where the phase is single valued. For 
values of M > a the function d{f3) is multivalued, allowing 
for the existence of many metastable static configurations 
at fixed u. Figure |SJ| shows one such multivalued S{P). 
Because of the metastability, the coherence can vary over 
some range. For a fixed u, the range in coherence results 
in a range of couplings fi. When u > a and (5(/3) is multi- 
valued, one chooses the (stable) branch of the (3{S) curve 
that is consistent with the particular metastable state 
one wishes to describe and also ensures that = 0, or 
equivalently that Eq. H10|l is satisfied. For simplicity and 
correspondence with "typical" sample preparation, we fo- 
cus on those metastable states accessed by adiabatically 
increasing u from zero*^ These correspond, for a given 
u, to the solid portions of the curve shown in Fig. ((S)). 
The details of the calculation for the scenario of adiabat- 
ically increasing u, which selects one branch, are given in 
Appendix A. ft is relatively straightforward to show that 
for a given u these are the states which have the largest 
coherence. This selection of largest-M states is consistent 
with our numerical calculations. Note that the form of 
the S{(3) curve and the discussion of multiple solutions 
is formally quite similar to parts of the calculation for 
the purely elastic case, though the physical motivation is 
rather different.^ 

The behavior of the coherence as a function of /i is 
shown in Fig. for four pinning potentials (for histories 
where the effective coupling u is adiabatically increased.) 
As anticipated on the basis of the perturbation theory, 
for a hard pinning force (curve (a) of Fig. 0) the coher- 
ence is a single-valued function of /i. The system exists 
in the zero-r IS phase for fi < jj,u- At ^„ there is a 
continuous transition to the CS phase, with r growing 
continuously from zero. For soft pinning forces (curves 
(c), cubic pinning force, and (d), sine pinning force, of 
Fig. with c < the coherence is a multivalued func- 
tion of ^. In this case the IS phase is stable up to /i„ when 
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FIG. 5: A sample plot of 5, the displacement of a degree of 
freedom from the minimum of the pinning potential, versus 
the pinning phase I3n{5) for branch numbers —2 < n < 2. 
The solid portions correspond to even n, while the dashed 
portions correspond to odd n. The global phase ip is chosen 
to be zero. Here, the effective interaction is large enough, 
u > a, that 5(P) is multivalued. The maximum magnitude of 
5 is denoted by 5max. 



/i is ramped up from below. At /i„ the coherence jumps 
discontinuously to the stable upper branch of the curve 
corresponding to the CS phase. When is ramped down 
from above the system remains in the CS phase down 
to the lower value /i^. For this class of pinning forces, the 
IS-CS transition is always hysteretic at F = 0. In the 
marginal case of a piecewise linear pinning force (curve 
(b)), r jumps discontinuously at the transition, but there 
is no hysteresis. 

The coherence curves shown in Fig.lHlcorrespond to the 
metastable states that would result through adiabatically 
increasing u. As mentioned earlier, for a given w, this 
is the state whose phases are as close as possible to the 
global phase V' = 0, and hence is the state with the largest 
coherence. Thus, the curves shown in Fig. |S1 are upper 
bounds on the coherence for each type of pinning force. 
In order to calculate the lowest possible coherence at each 
u, one must consider the metastable state whose phases 
are as far as possible from the global phase. To obtain 
this lower r(/.t) bound analytically is tedious, and we have 
done so only for the sawtooth linear case. This result for 
the lower bound is displayed in Fig. along with the 
upper bound, which, again, is the relevant state for the 
histories we consider here. 

In addition to determining the transitions between the 
IS and CS phases, the non-perturbative treatment at zero 
drive can also be used to determine if there is a critical 
value of /i, /XT, above which the depinning threshold van- 
ishes and the system is always sliding for all F > 0. 
We present an outline of the argument here and relegate 
the details of the calculation of ht to Appendix A. The 
threshold force can be thought of as the largest value of 
the driving force at which there still exists a stable static 
solution to the equation of motion. All such solutions 



10 




FIG. 6: The coherence of the static state at F = as a 
function of the coupling strength for four pinning forces: (a) 
hard (c > 0) cubic pinning force, with a = c = l/(7r + 7r^); (b) 
piecewise linear pinning force, with a — 1/n; (c) soft (c < 0) 
cubic pinning force, with a = l/vr and c — — I/tt"^; (d) sine 
pinning force whose maximum strength is 1/n. Also shown 
is the value fid where the coherence jumps from a finite value 
to zero upon decreasing jj,. 
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FIG. 7: Upper and lower bounds for r{u), plotted as the co- 
herence r{fi = u/r), corresponding to the maximal and min- 
imal coherence static metastable states. The pinning force is 
taken to be piecewise linear with a — 1/tt. A single static 
coherent solution is obtained only in the limit /i ^ oo. 



satisfy the static self- consistency condition. For incoher- 
ent static solutions, in which the domains are completely 
decoupled, this threshold force is simply the single par- 
ticle depinning force. For coherent static states the solu- 
tion 6{(3) is multivalued, but only those metastable states 
which satisfy the imaginary part of the self-consistency 
condition are acceptable solutions. Consider a system 
in which there are multiple metastable static solutions 
at zero drive. When an infinitesimal driving force is 
applied a correspondingly infinitesimal number of these 
states becomes unstable as they no longer satisfy the self- 
consistency condition. The system remains, however, 
pinned provided there still exist other accessible static 
metastable states. As the force is further increased, more 
static states become unstable, but the system does not 
depin until the "last" of the available static solutions, 
that is the one corresponding to the largest value of F 
for which a metastable static state exist, becomes unsta- 
ble. This value of F defines the depinning threshold. On 
the other hand, if there is a unique metastable static so- 
lution at zero drive, the system will depin immediately 
upon an infinitesimal increase of the driving force. When- 
ever there is a unique solution at = 0, the depinning 
force is therefore zero. As shown in Appendix B, for dis- 
continuous forces there are always a variety of metastable 
static states at zero drive for any finite value of /i (see 
also Fig. 01 J so that fir — oo. For continuous pinning 
forces, there is a finite coupling fix above which there is 
a single static state at zero drive and where the threshold 
force vanishes. This is for instance the case for the sinu- 
soidal pinning force, where the upper and lower bounds 
of r(^) (shown in Fig. 0) coincide and /ir = /^«- For a 
general continuous pinning force /xt is given by 



7r|r'(7r)| 



j;dd^l-iY{S)/Y'in)y 



(23) 



IV. STABILITY OF THE STATIC 
INCOHERENT PHASE AT NONZERO DRIVE 

We next consider static states in the presence of a fi- 
nite driving force, F ^ 0, starting with incoherent static 
solutions. We will use a perturbative treatment analo- 
gous to that of Sec. IIIII to analyze the limit of stability 
of the IS phase against varying fi and F. For finite F, 
the IS phase can become unstable to either the coher- 
ent static phase or the moving phase. The perturbative 
analysis described in this section allows us to establish 
whether the transition from the IS to CS phase at finite 
F is continuous or hysteretic, in much the same way as 
done in Sec. IIII Bl for F = 0. Again we find that the 
nature of the transition depends on the type of pinning 
potential, but the addition of a driving force changes the 
shape of the effective pinning force. This change can, 
in some cases, change a continuous IS^CS transition at 
F = to a hysteretic transition at finite F. The value 
of F above which the CS phase becomes unstable to a 



11 



moving state cannot be determined perturbatively and 
we defer its calculation to the next section. 

The perturbation theory described below is of course 
only valid for forces less than the single particle depinning 
force, Fsp. This force is the maximum value of 1^(2;) 
and is the driving force required to set in motion a single 
independent domain. It is hence the threshold force for 
an incoherent group of domains. 

We will study the stability of the incoherent phase to 
small changes in the coherence r. Taking the initial static 
phase to be incoherent, the effective coupling u = fir = 
and the static solution is obtained by simply balancing 
the pinning and driving forces. From Eq. 10} (with ^ = 
and h = 1) the noninteracting static solution is 9 — (3 — 
Y~^{F). It is convenient to choose the global phase to be 
non-zero, ip = —Y~^{F), and to work with the deviation 
6 — 9 — P + Y^^{F) from the incoherent static solution 
at a given F. The static solutions are then given by 



= F - usm{S + (3) + Y{S + So) , 



(24) 



where Sq ^ = —Y ^{F). For small u we can expand 
the pinning force in powers of S, 



YeffiS) 



Y{S + So)+F 
a{F)S + b{F)6'^ 



(25) 



c{F)S^ 



The effective pinning force Yeff{S) has precisely the same 
form as Y{6) for zero F, but the coefficients now depend 
on F through Sq = Y~^{F). These modified coefficients 
are given by 



a(F) = Y'iSo) , 
6(F) = Y"{So)/2 , 
£{F) = Y"'{6o)/6 



(26a) 
(26b) 
(26c) 



At non-zero drive the coefficient b{F) is always finite, 
reflecting the fact that the external drive makes the pin- 
ning force asymmetric about ^o- The equation for (5(/3) 
is then formally identical to that for 5{(3) in the F = 
case, with Y{S) Yeff{S), 



= -ufim{6 + P)+Yeff{6) . 



(27) 



Similarly, the self consistency conditions can be expressed 
in terms of S as 

f d/3cos(5 + /3)=/(u,F), (28) 
where r is now a function of both u and F, and 



= / d/3sm{S + l3) . 



(29) 



We can now use the results obtained in the zero drive 
perturbation theory. The value of fi at which the IS phase 
becomes unstable is given by 



Ai„(F) = 2&{F) , 



(30) 



and will now in general depend on F. Conversely, we can 
define a critical line Fu{n) as the solution of /i = 2a (F^). 

For drives sufficiently small that the system remains 
pinned at the instability line, the form of the onset of 
coherence near fJ.uiF) can be determined by looking for 
a solution to Eq. (|28|l in the form of a power series. 



(31) 



f{u,F) :^ri{F)u + r2{F)u^ +r3{F)u^ + 



As usual in such calculations, we expect the nature of the 
instability to depend on the signs of the coefficients. The 
coefficients ri{F), r2{F) and r^{F) are given by Eq. (|18c() 
with a, &, c replaced by a{F), b{F), c{F), giving 



ri{F) 
r2{F) 



1 



2a{F) 




3 ( h{F)i{F)~2h{Ff 
8 I h{Ff 



(32a) 
(32b) 

(32c) 



Thus, the form of r(//, F) near fJ.u{F) is: 



r(u F) = 1 (ji^MEL 



1/2 



(33) 



As for the case of F = 0, the behavior is controlled by the 
sign of the coefficient »'3(F) of the cubic term in Eq. H31|l . 
If r3(F) > the coherence grows as ~ (/i — /x„(F))^/^ 
with increasing /x, indicating that the r versus fj, curve is 
continuous. Conversely, if rsiF) < the coherence grows 
with decreasing /i as (/ii,(F) — /i)^/^, and the r versus fi 
curve is hysteretic. One important complication is that 
for finite F the coefficient r^ can change sign as a function 
of F for a given pinning force. As a result the transition 
between coherent and incoherent static states can change 
from continuous to hysteretic above a characteristic force 
Fh defined by the solution of r3(Fh) = 0. 

We now specifically apply these general results to the 
three classes of pinning forces (linear, hard, and soft.) 
Again, these are of the general form 



Yix) 



-ax — cx 



'TT < X < TT 



(34) 



with a > 0. The three classes have c zero, positive and 
negative, respectively. 



A. Piecewise Linear Pinning Force (c = 0) 

For the piecewise linear pinning force of Fig. ^e), 
where Y{S) is given by Eq. H34|l with c = 0, we sim- 
ply have a(F) = a and 6(F) = c(F) — 0. In this case 
i^u{F) = fiu{0), independent of F. In fact we will show 
in Section lV Al that the coherence r(/i) of the entire static 
state is independent of F for all values of /z, whenever the 
system is pinned. The IS phase is stable for /i < /i„ = 2a 
and F < Fgo = air. 



This region is shown in Fig. |S1 
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FIG. 8: The region of stability of the IS phase for a piecewise 
Unear pinning force. The single particle depinning force Fsp 
and the coupling strength Hu for instability to the coherent 
or moving states are also indicated. 




FIG. 9: A plot of the region of stability of the incoherent 
static phase for a hard cubic pinning force. The nature of the 
instability along the Fu{i-i) curve is indicated by the thickness 
of the bounding curve on the right. For jj. > /^h (F > -fh), 
the transition is hysteretic, while for smaller couplings (or 
small, fixed driving force for varying couplings) the transition 
is continuous. 



B. Hard Cubic Pinning Force (c > 0) 



In Fig.|21we show the region of stability of the IS phase 
for the hard pinning force of Fig.^d). In this example, 
the maximum pinning force from Eq. (|34|l gives the single 
particle depinning threshold as Fgp = ott + cn^. When 
the coupling fi is ramped up adiabatically with constant 
F < Fsp, the IS state becomes unstable at a value /x„ 
given by (see Eq. (jSnjl), 

^luiF) = 2[a - 3c6liF)] . (35) 

For the hard cubic potential this result can be inverted 
analytically to obtain the boundary F„ (/i) of the IS state 
shown in Fig. |51 with the result 



Fu{^^) 



6 



6|c| 



(36) 



The maximum value of /i for which the IS state is stable 
is /X*, where fi* is found by the intersection of the IS 
depinning curve and the Fu(^jl*) curve. Its value is 



6cn 



(37) 



Note that if the system is prepared in the IS state at > 
/Lt„, then a transition to a coherent state can be achieved 
by decreasing F. This is because decreasing F allows 
the domains to relax back toward the minima of their 
pinning potentials, where the pinning force (determined 
by the curvature of the potential) is smaller and hence 
the coherence can increase. 

In the case of the hard cubic pinning force the coeffi- 
cient r3{F) can change sign as a function of F. For small 
F, r3{F) > and the transition from the IS to a coher- 
ent static phase is continuous. Above a critical value Fh 
defined by r3{Fh) = the transition becomes hysteretic. 
The force Fh is given by 



16 

15372 



(38) 



and is small compared with Fgp for the potential shapes 
and parameters we have considered. For a = c ^ 1/{tt + 
TT^), we find Fsp = 1, Fh « 0.008 and « 1.2m„. 



C. Soft Cubic Pinning Force (c < 0) 

Soft cubic pinning forces given by Eq. (|34|l with c neg- 
ative, can be divided into two classes: (i) forces that are 
monotonic functions of the phase within each period, as 
plotted in Fig.^a), and (ii) those that reach their maxi- 
mum (minimum) within a given period and turn over, as 
plotted in Figs.^b) andQIc). Holding /i constant, the 
incoherent static state becomes unstable upon increasing 
F to Fuifi), with 



Fuip) 



2/i„ + n ffl 



6 



6 c 



(39) 



unless the single particle depinning force is first reached. 
For pinning forces in class (i) the value /i* where 
Fu(/i*) — Fsp is positive and the region of stabihty of 
the incoherent static state is of the type shown schemat- 
ically in Fig. llOf a'! . For pinning forces in class (ii) (for 
pinning forces with only cubic terms, this class is given by 
c| > 1/{3tt^)), it can be shown that fi* = 0. The single 
particle depinning transition is always preempted. Here, 
the region of stability of the incoherent state is deter- 
mined by Fu{n) for all values of /i, as shown in Fig. lior bb 

For a soft cubic pinning force r^{F) is negative for all 
F and the transition from the incoherent to a coherent 
state is always hysteretic. 



13 




FIG. 10: Sketches of the region of stabihty of the incoherent 
static phase for a soft cubic pinning force. Figure (a) corre- 
sponds to a pinning force of type (a) that does not turn over 
(is monotonic) in each repeated interval. Figure (b) corre- 
sponds to a pinning force of type (b) that are non-monotonic 
in each period. 



V. NONEQUILIBRIUM PHASE DIAGRAMS IN 
THE fi-F PLANE. 

In this section we present the nonequilibrium phase 
diagrams in the fi-F plane for the various pinning forces 
introduced in Fig.^ The phase diagrams are based upon 
both analytical results and numerical computations. The 
analytical bounds on the stability of the static phases are 
based on the previous sections' results for the incoher- 
ent static phase and calculations for the coherent static 
phase whose details are presented in the appendices. Nu- 
merical integration of the equations of motion is used 
to determine the boundaries of the moving phases: by 
starting from the moving phase and decreasing F or fi, 
the repinning curves can be found. Of special interest 
is the nature of the depinning transition obtained when 
the applied force F is varied at constant fi. The curves of 
mean velocity as a function of driving force correspond to 
the IV characteristics of physical systems, such as CDWs 
and vortex lattices. Our focus is on classifying models or 
parameter ranges for which the depinning transition is 
continuous or hysteretic. In general, for each of the pin- 
ning forces we consider, the depinning transition appears 
to be continuous with a unique depinning threshold at 
large ^, where the system is more rigid. In contrast, the 
velocity-force curves generally exhibit macroscopic hys- 
teresis at small values of /i, where the system is more 
likely to display plastic effects. 



A. Piecewise Linear Pinning Force 

In Sec. IIV Al perturbation theory was employed to 
study the transition between incoherent and coherent 
static phases for the piecewise linear pinning force. It 
was found that when the coupling strength fi is changed 
at fixed F within the pinned region of the phase diagram 



this transition is always discontinuous, although not hys- 
teretic. Furthermore, the critical value of fj, where the 
transition takes places appears to be independent of the 
driving force. Here we show that this remains true in 
a complete calculation. We also calculate the depinning 
threshold exactly by determining the limit of stability of 
the static phases. For the piecewise linear pinning force 
(i.e., Eq. with c = 0), the force balance equation in 
the static state is 

= F-u sm{S + f3)-a5, (40) 

where — tt < 6 < n and we have chosen -0 = 0. Letting 
5 = S + F/a and (3 — a ~ F/a, Eq. (|in|) can be written 
as 

= -usm{S + a) - a5 , (41) 

with — TT — F/a < 5 < IT — F/a. It is apparent from 
Eq. H41I) that (5 is a function only of a and u and does 
not depend on F explicitly. The real part of the self- 
consistency condition that determines the coherence r 
becomes 

1 r 

r=— dacos{S + a), (42) 

and clearly r{u, F) = r{u, F — 0). Thus, the coherence of 
the static state is independent of F. The line separating 
the incoherent and coherent static phases is a vertical line 
a,t ^ = = 2a in the fj, — F plane, as shown in Fig. 1111 
The IS-CS transition is discontinuous and non-hysteretic 
at all values F where the static phases are stable. When 
the force is ramped up adiabatically at fixed fi < Hu from 
the IS phase where r = 0, the system depins at the single 
particle depinning force Fsp = air. For ii > fj.u the system 
is in the CS phase, where the coherence is non-zero and S 
is a multivalued function of a. As discussed in Sec. lIII CJl 
there are many static metastable states available to the 
system for a fixed value of u. We relabel the metastable 
states and denote each state by a (Si (a, u) which is a single 
valued, but generally discontinuous, function of a. Each 
6i must satisfy the imaginary part of the self-consistency 
condition which using Eq. (|41|l can be rewritten as 

= /" daS,{u,a) . (43) 

This implies that the acceptable Si's are odd functions of 
a. In addition, each static metastable solution must lie 
within the upper and lower bounds, Su{F) = n — F/a and 
Si{F) = — TT — F/a. As F is increased, the value of the 
upper bound decreases, reducing the number of allowed 
(5i's, until at = F^{ii) only one solution remains. This 

special state, S{a,u), is equivalent to the one that would 
be obtained through adiabatically increasing u. The as- 
sociated r(/i) curve is shown in Fig. |Blfor a = I/tt. The 
value of F|(u)S& is given by Smaxiu) = tt — F^/a. For 
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FIG. 11: Phase diagram for the piecewise hnear pinning 
force, Y(x) = —x/n (see Fig. Q^))- The hghtly shaded 
portion is the coexistence region of the IS and CM phase 
(/i < yLu) and the smaller, darkly shaded region, is where the 
CS and CM phases coexist (p^ < fJ, < fJ-e)- The depinning 



lines F| = F^p 



and F| have been obtained analytically and 



confirmed by numerics. The boundary where the system 
repins was obtained numerically. The point {^^,Fe) marks 
where the static-moving transition changes from hysteretic to 
continuous. The boundary between the IS and CS phases is 
F independent and lies at /i = 



u < an/ 2 we find from Eq. H41|l 5„ 
gives 



.(u) — u/a which 



u < a7r/2 



(44) 



For u > a7r/2 Smax is defined implicitly by aS„ 
usm{Smax) and is given by 



u sm 



[irFfifi) -u/a] , u>an/2. (45) 



It is then possible to calculate F^(p) using the expression 
for r{u) given in Eq. ljA12|) . The resulting phase diagram 
is shown in Fig. ITTI 

For < the static phase is incoherent and the de- 
pinning transition is hysteretic in both v and r, as shown 
in the top two frames of Fig. ^| The system depins 
at Fgp when the drive is ramped up adiabatically from 
the static phase, but repins at the lower force F^ when 
the force is ramped back down from the sliding state. 
The line F^ has been obtained by numerical simulation 
of the mean field model. The numerics have also revealed 
that a small region of hysteresis persists for ^ > al- 
though the static phase is coherent here. The behavior 
of V and r in this region is shown in the two mid frames 
of Fig. El Finally, for fj. > fie (where /ig is the value of 
the coupling above which the static-moving transition is 
elastic in nature) the depinning is continuous, as shown 
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FIG. 12: Mean velocity v and coherence r as functions of the 
driving force F for the piecewise linear pinning force. The 
curves are obtained numerically by first ramping F from zero 
to a value well within the sliding state (F = 1.2), and then 
decreasing F back down to zero, while holding fi constant. 
The top frames show the behavior for /i = 0.25, where the 
initial static state is incoherent: this state starts sliding at 
the single-particle depinning force F^p — 1 and repins at a 
lower force F « 0.88. The mid frames display the results for 
an initially coherent static state (/x = 0.64 > fj.u), which still 
displays hysteresis, both in v and r. The bottom frames are 
for fj, = 1.0, which has an initial coherent state and undergoes 
continuous depinning. 



in the bottom frames of Fig. E| The values of /ig and F^ 
are defined via F^{fie) = Fi{fie) — Ff.. Finally, the de- 
pinning threshold F^ is nonzero and finite for all /i, i.e., 
fiT ~ oo. This is a general property of discontinuous pin- 
ning forces, to be contrasted with the behavior observed 
for continuous pinning forces, such as the sinusoidal one 
studied by Strogatz and coUaboratorsi^^ 

Before closing this subsection, we must address the 
possibility an of incoherent moving (IM) phase. Strogatz 
and collaborators^^ found that an IM phase is always 
unstable for a sinusoidal pinning potential. It can be 
shown that this remains true for other continuous pin- 
ning forces. The situation is less clear for discontinuous 
pinning forces. In Appendix D we present the details of a 
short time {t = Q) stability analysis for the IM phase for 
any Y{x). This analysis will tell us something about the 
long time, steady state limit, provided r{t) is a monotonic 
function of time. This analysis predicts a range of sta- 
bility for the IM phase for discontinuous pinning forces, 
provided the jump discontinuity at a; = tt is taken into ac- 
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FIG. 13: Phase diagram in the couphng-drive (i^-F) plane 
for a hard cubic pinning force of the type shown in Fig. |^) . 
The form of the pinning force Y{x) is given by Eq. (1341 . 
with a = c = l/(7r + 7r^). The regions of IS-CM, CS-CM 
and IS-CS-CM coexistence are shown in light, medium and 
dark gray, respectively. The incoherent and coherent depin- 
ning lines are denoted by and Ff respectively. The re- 
pinning line is denoted by F^. The coherent depinning line 
and the repinning line join at (/ie, Fe). Beyond this point the 
static- moving transition is continuous. The curves Hu{F) and 
^id{F) are the values of the coupling at which the static sys- 
tem makes the transition to and from finite coherence states, 
respectively. There curves join at (^h,-fh) where the IS-CS 
transition becomes continuous. 

count when preparing the system. However, simulations 
show that r(t) is in general not monotonic and that the 
strength of the perturbation needs to be decreased with 
system size in order to observe the IM phase, suggesting 
that the perturbative short-time analysis is simply not 
valid in this case. Finally, if a narrow distribution of pin- 
ning strengths h is introduced, we find numerically the 
IM phase to be unstable. Given these numerical find- 
ings, we believe that the IM phase is generally unstable 
in mean field theory. 



B. Hard Cubic Pinning Force 

The phase diagram for a hard cubic pinning force, 
given by Eq. I|34() with c > (see Fig. H^d)) is shown 
in Fig.inifor a = c= 1/{tt + t:'^). 

Though the general topology is similar to that of the 
phase diagram for the piecewise linear force, the history 
dependence is significantly more complicated. A first 
difference is that the transition between the IS and CS 
phases is now continuous for F < Fh, with given by 
Eq. (|38|l . and hysteretic for F > Fh. For the parameter 
values displayed in Fig. E| the value of Fh is very small. 



but still finite. A second new feature of the phase dia- 
gram is the presence of a small region (darkest gray in 
Fig. 1131) where all three phases coexist. 

The strong history dependence is manifested in the 
macroscopic response and includes reentrant behavior for 
fixed fi or F histories. The mean velocity and coherence 
are plotted as a function of (increasing, then decreasing) 
driving force for a few typical values of /i in Fig. 1141 The 
pinning force is given by Y(x) = — (a; + x'^)/(7r + 7r''). The 
top frames show a simple hysteretic depinning transition 
for a system prepared in the incoherent static state at 
F — 0, similar to that seen for a linear pinning force. 
The middle row of frames display the more complicated 
history that results when the system is prepared in a 
coherent static state at F = 0, with ^ = 0.5. The 
velocity shows a single hysteresis loop, but the plot of 
coherence r shows first a decrease and then a jump to 
the incoherent state as the force is increased, followed 
by a jump back to a finite value when bulk depinning 
takes place. In this case both the regions of IS-CS and 
IS-CM coexistence are crossed when F is ramped up. 
The IS-CS transition occurs as the phases are pushed 
away from their zero-force minima to regions of the pin- 
ning potential with higher curvature, which makes the 
coherent state unstable. Upon decreasing the force, both 
the coherence and the velocity jump back to zero, then 
the coherence increases again as the force is decreased. 
The jumps in coherence when F is ramped down occur 
at values of F different from those where the coherence 
jumps during the ramp up. For rather specific values of /i, 
even more baroque histories can be found by crossing the 
three-phase coexistence regions. An example is shown in 
the last row of frames in Fig.^l where /i — 0.76. Here, 
the sequence is CS— s-IS^CM^CS, which skips the IS 
phase on decreasing F. Note that the velocity vs. drive 
force curve is relatively unremarkable, showing simple 
hysteresis in this case. The coherence history is more 
complicated. 

Another interesting feature of the phase diagram for 
the hard cubic pinning force is that at constant fi, a por- 
tion of the moving phase lies between the incoherent and 
coherent static phases. This suggests the possibility of 
re-entrance in the depinning transition for fi > fXe- It is 
not, however, straightforward to prepare the system in 
the lightly shaded portion of the phase diagram where IS 
and CM phases coexist and n > fie The static solution 
must either be created "by hand" at that location {fi, F) 
in phase space or the system can be prepared in the IS 
phase at a lower value of fi and the coupling can then 
be ramped up to the relevant value fi > fi^. Both the 
difficulty of preparing the system in the re-entrant state 
and the re-entrance for a specially prepared state are dis- 
played in Fig. El Here both sets of curves correspond to 
the same value of the couphng strength, fi = 1.25. In the 
top pair of curves the system is prepared in the coherent 
state at F = 0. As the force is ramped up adiabati- 
cally, the system depins continuously at F^ , where both 
velocity and coherence change smoothly, with r rapidly 
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FIG. 14: Mean velocity and coherence versus force for the 
hard potential and various values of /i. Solid lines are used to 
display the response obtained when F is ramped up from zero, 
while dashed lines show the jumps in v and r when ramping F 
back down. The top frames show the hysteretic depinning of 
a system prepared in the IS phase. For /i — 0.5 (mid frames) 
the system is initially in a coherent (r 7^ 0) static state at 
F = 0. As the force is ramped up, the system first crosses 
the boundary from the CS to the IS phase, where r jumps 
discontinuously from its initial finite value to zero, while the 
system remains pinned {v = 0). At a higher force the system 
depins by crossing the boundary from the IS to the CM phase 
and r jumps from zero to a large finite value. The subsequent 
ramping down of the field goes through this sequence of phases 
in reverse order, but the jumps occur at distinct values of F. 
The bottom frames describes the complex response that takes 
place along a path that crosses the dark region of three-phase 
coexistence. See the text for further description. 



approaching its limiting value, r — 1. The coexistence 
region is never accessed in this case. In the bottom set 
of figures, the system is prepared in an incoherent static 
state at finite F, deep inside the coexistence region. The 
system is then observed to depin as the force is ramped 
down at constant /x across the boundary between the co- 
existence region and the CM phase. Simultaneously, the 
coherence jumps from zero to a large finite value. Upon 
further ramping down F, the system repins again contin- 
uously at F^. 



C. Soft Cubic Pinning Force 

We distinguish three types of soft cubic pinning forces 
^iven by Eq. H34|l with c < 0. These pinning forces and 
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FIG. 15: Both sets of figures show the behavior of velocity 
and coherence for /i = 1.25, but for different initial states. 
The top frames are obtained by preparing the system in a 
coherent state at /i = 1.25 and -F = 0, and ramping F up to 
a value above F^, and then back down to zero. In this case 
the depinning is continuous. The bottom frames are obtained 
by preparing the system in an incoherent state a.t ^ = 1.25 
and F = 0.9, inside the lightly shaded area of coexistence of 
CS and CM phases, and then ramping the force down to zero. 
Note the depinning upon decreasing force in this case and the 
subsequent repinning. 



corresponding potentials are shown in Fig. ^ (a) forces 
that are monotonic over the entire period and do not turn 
over in the interval [— tt, tt]; (b) forces that are nonmono- 
tonic over the period and do turn over in the interval 
[— TT, tt], but are discontinuous; and (c) continuous forces, 
which are obviously nonmonotonic. The phase diagrams 
for these potentials exhibit qualitative differences as com- 
pared to those discussed so far. Specifically, the CS re- 
gion at non-zero F may or may not extend to /x = 00 
and may not even exist. For most potentials, however, 
we do find a non-trivial coherent static phase. The only 
exception is the case of a sinusoidal pinning force studied 
previously by Strogatz and coUaboratorsjSSi where the CS 
state is unstable. 

For monotonic pinning forces, (a), the boundaries -F| 
and i^sp intersect at a finite positive value /i, of /x, given 
by Eq. (|37|l (see Sec. IIV Cl for a full discussion). This re- 
sults in a portion of the depinning boundary being hor- 
izontal on the ^-F plane, as F^ — Fgp for fj, < fj,^, as 
shown in Fig. 1161 In contrast, if the pinning force is non- 
monotonic, (b) or (c), and reaches its maximum within 
the period, then /i* = and the phase boundary has no 
horizontal portion. This behavior is shown in Fig. [3 for 
a nonmonotonic, but discontinuous pinning force. 

The results for pinning forces of type (c), that are con- 
tinuous (and therefore must be non-monotonic) have two 
important features: fi* = and fiT is finite. These fea- 
tures imply, respectively, that there is no horizontal por- 
tion to the CS depinning curve and that the system slides 
at arbitrarily small force whenever the coupling is large. 
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FIG. 16: Phase diagram in the coupling-drive (fi-F) plane 
for a soft monotonic cubic pinning force of the type shown in 
Fig. ^a). The pinning force Y{x) is given by Eq. 134^ with 
a = 6/(57r) and c = -I/Stt^ The regions of IS-CM, CS-CM 
and IS-CS-CM phase coexistence are shown in hght, medium 
and dark gray respectively. The lines and are the forces 
at which the system depins upon increasing the drive from the 
incoherent and coherent static states respectively. The line 
Fi is the force at which a moving system stops upon lowering 
the drive. The coherent depinning line and the repinning line 
join up at (/ie,Je) and the static-moving transition becomes 
continuous. The curves Hu{F) and ^d(F) are the values of 
the coupling at which the static system makes the transition 
to and from finite coherence states respectively. 



i.e., when /i > ^lt- The typical phase diagram for a pin- 
ning force of this type is shown in Fig. El Figure ITSl 
shows sample v{F) and r{F) plots for this case. 

At finite drive the CS phase does not extend beyond 
/i = 1.84. For values of the coupling between /z^ and 
the CS phase exists at finite drive, albeit only for very 
small values of F < F^{)jl) <g; 1. This small region of the 
phase diagram in Fig. 1171 is magnified and shown in the 
inset. It is interesting to compare the these results with 
those obtained by Strogatz and collaborators^^ for an- 
other continuous pinning force, namely Y[x) — — sin(a;). 
The corresponding phase diagram is shown in Fig. 1191 In 
this case fiu = fJ-T = 2 and, more significantly, F^{fi) — 0. 
This means that the CS phase never exists at finite F. 
Thus, it seems that sinusoidal pinning forces are a spe- 
cial class of more general continuous pinning forces in 
that they never allow the possibility of a CS phase at fi- 
nite drive. This difference, while important qualitatively, 
may not be quantitatively significant given that F^(fj,) is 
always very small. 

Finally, for any continuous pinning force, the IM phase 
is not stable even in the short time analysis. (See Ap- 
pendix D.) This result is consistent with the findings of 



FIG. 17: Phase diagram in the coupling-drive (i-i-F) plane for 
a soft cubic pinning force of the type shown in Fig. H^c) . The 
pinning force Y{x) is given by Eq. (O with a = 3^3/ (27r) 
and c = —3^/3/{2tv^). This choice of parameters gives a non- 
monotonic and continuous pinning force: the results are to 
be compared with Y{x) = — sin(a::), another non-monotonic 
and continuous force. The region of IS-CM coexistence is 
shown in light gray, while the IS-CS coexistence is shown in 
medium gray. The lines F^ and F^ are the forces at which the 
system depins upon increasing the drive from the incoherent 
and coherent static states respectively. The line Fi is the 
force at which a moving system stops upon lowering the drive. 
The CS region is very small for these values of parameters, 
corresponding to /it ~ 1.85 and fJ,d{0) ~ 1.44, and it has been 
magnified in the inset. The CS phase does not exist at finite 
F for coupling larger than ^t- Shown within this inset is the 
point {fj.e, Fe) where the Ff and Fi lines join and the static- 
moving transition ceases to be hysteretic. The curves fJ,u{F) 
(only visible within the inset) and ^id{F) are the values of the 
coupling at which the static system makes the transition to 
and from finite coherence states respectively. 



Strogatz and collaborators^^ as well as our simulations. 



VI. AVERAGING OVER DISORDER 

In this section we discuss the role of the shape of the 
distribution p(h) of pinning strengths on determining the 
nonequilibrium phase diagram. In the previous sections 
we restricted ourselves to an infinitely sharp distribution, 
p{h) — 5{h — 1). This choice is appropriate for systems 
with strong pinning and allows for a direct comparison 
with the results of Strogatz and collaborators^^. It is 
easy to show that the nonequilibrium phase diagram of 
the driven system retains the same qualitative structure 
for any distribution that is sharply peaked around a fi- 
nite value of the pinning strength and vanishes below a 
finite ho > 0. A broad distribution of pinning strength 
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FIG. 18: Mean velocity and coherence, obtained from nu- 
merical calculations, for a continuous cubic pinning potential 
and parameter values given in Fig. 1171 The top frames record 
the hysteretic response of a system prepared in the incoherent 
state at fi — 0.8 and F = 0, while the bottom frames show 
the continuous _F = depinning of system prepared in the 
coherent state a.t fi — 1.67. 



may, however, qualitatively alter the mean field physics. 
Broad distributions p{h) are of interest to model phys- 
ical systems with weak pinning. Furthermore, a broad 
distributions of pinning strengths yields variations of the 
local stresses in the mean field theory and may give us 
some insight into the behavior of the system in finite di- 
mensions. 

We consider a distribution of pinning strengths p{h) 
that vanishes below a minimum pinning strength ho > 
0. As will become apparent below, it is important to 
distinguish three class of distributions: 

1. distributions that vanish below a finite pinning 
strength, i.e., p{h) = for h < ho, with Hq > 0; 

2. distributions with no finite lower bound of the pin- 
ning strength, but zero weight at h — 0, i.e., ho — 0, 
and p{0) = 0; 

3. distributions with no finite lower bound of the pin- 
ning strength, and finite weight at h = 0, i.e., 
ho = 0, but p{0) > 0. 

The nonequilibrium phase diagram depends qualitatively 
on whether or not the lower bound Hq is finite. If the 
distribution of pinning strengths p{h) vanishes below a 
minimum pinning strength Hq > 0, the single particle 
depinning threshold Fgp remains finite and the system 
exists in an IS phase for F < Fgp. When ho = 0, the 
single particle depinning threshold vanishes and the IS 
state can only be stable at = 0. 

If the IS phase exists, its stability can be analyzed for 
an arbitrary distribution p{h) by the perturbation theory 
described in Sec. IV. For arbitrary h, the static force 



FIG. 19: Phase diagram for a sinusoidal pinning force 
Y{x) — — sin(a;). The IS-CM coexistence region is shaded 
gray. The F — region in which the system can only exist 
in the CS phase is denoted by a series of x's. The region of 
IS-CS coexistence is denoted by medium on-axis gray shad- 
ing. The IS^CS and IS^CS phase boundaries are the points 
{p = Pu = 2, F = 0) and {p — pd ~ 1.49, F — Q) respectively. 
The line F^ is the forces at which the system depins from the 
incoherent state. The line Fy is the force at which a moving 
system stops upon lowering the drive. 



balance equation has the form 

= F - M sin(6l - 1/^) + hY{e - f3) , 



(46) 



with the self-consistency condition given by Eq. I^. 
Clearly this equation is identical to the equation stud- 
ied in Sec. IV for h — 1, provided we rescale both the 
driving force F and the coupling strength u by the pin- 
ning strength h. We can then carry out the perturbation 
theory described in Sec. IV as a perturbation theory 
in powers of u/h, provided of course u << ho- This 
shows that the perturbation theory breaks down when 
ho 0. Furthermore we must require F < Fgp, which is 
a necessary condition for the existence of the IS phase. 
Proceeding precisely as in Sec. IV and using the same 
notation, we obtain an expression for the coherence r as 
a power series in u/h, given by 



r = fiu,F) = 



dhpih)[ ry{F/h)- + r2 
h 



with 



n{F/h) 

r2{F/h) 
rsiF/h) 



2a{F/h) ' 
0, 

3a(f//i)g(f) -2[6(F/fe)]' 
8 h{F/hf 



(47) 

(48a) 
(48b) 

(48c) 
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and 



~a{F/h) = y'(5o) , 
h{F/h) = Y"{6o)/2 , 
c{F/h) - y"'(Jo)/6 . 



(49a) 
(49b) 
(49c) 



The boundary of stability of the IS phase, /i„(i^), is 
obtained hke before by solving the implicit equation 
r(^, F) = f{u = lJ.r,F) with f{u,F) given by Eq. (|?7|l . 
with the result 



1 



n -1 



ha{F/h) 



(50) 



If the distribution p{h) vanishes below a finite minimum 
pinning force ho > 0, then fXu remains finite and there is 
a range of /i and F where the IS phase is stable. Con- 
versely, if ho 0, the integral in Eq. (|50|l may diverge, 
yielding fi^ = 0. Below we will treat in detail the case 
of a piecewise linear pinning force, with Y(x) — ~ax. In 
this case Eq. ((301) reduces to 



2a 



dh 



p{h) 



(51) 



For concreteness, we consider a distribution of the form 

h> ho , 



p{h) - {h 
p{h) = , 



h < ho , 



(52) 



with ho > and a > 0. This form encompasses the three 
classes of distribution functions introduced at the begin- 
ning of the section. We can then obtain the boundary of 
the IS phase for a piecewise linear pinning force by eval- 
uating the integral on the right hand side of Eq. (|51|) . 
For distributions of the first class, corresponding here to 
ho > and a = 0, we find that Hu is finite at finite 
F and it is given by pu = 2ae^°Ei{ho), where Ei{x) is 
the exponential integral. For this type of distribution it 
can be shown that the nonequilibrium phase diagram re- 
mains qualitatively similar to the one obtained for the 
sharply pinned distribution, p{h) = 5{h — 1), even for all 
types of pinning forces studied in Sec. V. When ho ^ Q 
the perturbation theory breaks down and the existence 
of a finite value of even at = 0, depends on the 
form of p{h) for /i ^ 0. For distributions of the second 
class, with ho = 0, but p{Q) — 0, it can be shown that 
Pu is finite at F = 0, but vanishes at all finite F. In 
this case there is an IS-CS transition at F = 0, which is 
a remnant of the transition seen at finite F for the case 
of an infinitely sharp pinning strength distribution. For 
instance, for p{h) = he^^{a — 1), there is an IS-CS tran- 
sition at F = and pu — 0.27. Finally, for distributions 
in the third class, with p{Q) > 0, it can be shown that 
Pu vanishes as l/ln(l//io) when ho 0. For such distri- 
butions, there is no IS phase even at F = 0. The phase 
diagrams for this class of distributions of pinning strength 
are qualitatively different from those presented in Sec. V 
for all pinning forces. An example is shown in Fig. 1201 for 
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FIG. 20: Phase diagram in the p — F plane for a piecewise 
linear pinning force, with a — 1/tt and c = 0, and p{h) = 
e"''. The depinning curve has been obtained numerically for 
a system with A'' = 1024 and a ramp rate of dF/dt = 10"". 



the piecewise linear pinning force and p{h) — . This 
phase diagram has been obtained numerically. In the 
limit of large system sizes and adiabatically slow ramp 
rates dF/dt, no IS phase is observed even at F = 0. 
The small region of hysteresis in the transition between 
the CS and CM phases is also washed out by the disor- 
der averaging. The depinning curve F^ displays a broad 
maximum at a finite p and vanishes as p ^ oo. 

In general, the numerical simulations show that a 
broad distribution of pinning strengths with vanishing 
ho always washes out the IS phase and any hysteresis of 
the depinning transition. Whether this behavior persists 
in finite dimensions remains an open question. 



VII. DISCUSSION 

In this paper we have used a combination of analytical 
and numerical techniques to study the nonequilibrium 
mean field phase diagram of a model of an extended sys- 
tems with phase slips driven through disorder. For uni- 
form pinning, we generically find two stable static phases 
and a single moving phase. Both incoherent (IS) and co- 
herent static (CS) phases are possible, as well as regions 
where the two phases coexist. The moving phase, in con- 
trast, is always coherent (CM) in mean field theory. (An 
incoherent moving phase can be prepared by using special 
initial conditions, but does not appear to be stable.) Co- 
existence of two, or even three, of these phases can occur 
depending on the system preparation; this coexistence 
results in hysteretic transitions. Such a variety of phases 
was not found for the case of a sinusoidal pinning force 
analyzed earlier where only the IS and CM phases were 
found. While a discontinuity in the pinning force is not 
required for the existence of the new CS phase at large 
values of the coupling constant p, a jump discontinuity 
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FIG. 21: Phase diagram, redrawn in the disorder-drive plane, 
for a discontinuous soft cubic pinning force of the type shown 
in Fig. lb and pijn) — 5{h — 1). The disorder h and drive F 
are normahzed by the strength of the phase-slip interaction, 
fj,. The parameter values and the symbols are the same as in 
Fig. 2. 



in the pinning force does increase the range of F and /i 
over which the CS phase is observed. This is because 
discontinuity in the pinning force makes it more difficult 
for the system to depin, so that the static pinned phases 
can exist up to large coupling strengths, where the sys- 
tem is forced to acquire long range coherence. Once the 
system has become coherent, and therefore more rigid, 
the depinning threshold decreases with increasing ^, but 
remains finite for all finite values of the coupling strength 
and only vanishes for fi oo. For a continuous pinning 
forces, on the other hand, the depinning threshold van- 
ishes above a finite value of fi. 

In order to make some contact with particle simu- 
lations and with experiments, it is useful to discuss 
the mean field phase diagram in terms of the disorder 
strength h and the driving force F, rather than in the 
{n, F) plane as done so far. In most particle simulations 
it is the strength of the disorder that is most easily var- 
ied rather than the strength of the coupling. Disorder 
is also a crucial control parameter in many experimental 
systems. For instance, varying the applied magnetic field 
in current-driven vortex lattices has the effect of varying 
the strength of the disorder. At high fields the vortex 
lattice becomes softer and can better adjust to disorder. 
Increasing the magnetic field therefore corresponds to an 
effective increase of the disorder strength. Fig. [2] shows 
the mean field phase diagram in the {h, F) plane for the 
discontinuous soft cubic pinning force shown in Fig. lb. 
The corresponding phase diagram in the (/i, F) plane was 
shown in Fig. 2. 

When the disorder is weak relative to the strength of 
the coupling fi the static phase is coherent. At strong 



disorder the static phase is incoherent. The transition 
between the coherent and incoherent static phases at 
fixed /i is hysteretic with a region of coexistence of the 
two phases. At weak disorder there is a continuous 
"elastic-like" depinning transition from the CS to the 
CM phase. At large disorder the static phase is inco- 
herent and degrees of freedom depin independently at 
the single particle depinning threshold, F^. The moving 
system immediately acquires long-range correlations, be- 
coming much stiffer and harder to pin. As a result, when 
the force is ramped down the CM state repins at the 
lower force F^ . The qualitative features of this phase di- 
agram are remarkably similar to those obtained by Olson 
and coUaboratorsSS in a numerical simulation of a model 
of a current-driven layered superconductors, with mag- 
netically interacting pancake vortices. At weak disorder 
these authors find that the layers are coupled and the 
system forms a coherent three-dimensional static phase, 
with long-range correlations along the direction normal 
to the layers, which depins continuously. At strong disor- 
der the static state consists of decoupled two-dimensional 
layers. When the driving force is ramped up from this 
incoherent static state, the layers depin independently at 
the single-layer depinning threshold and the transition is 
hysteretic. One difference between our mean field model 
and the numerical model studied by Olson et al. is the 
absence, in our model, of an incoherent moving phase. In 
the layered superconductor at strong disorder the layers 
remain decoupled upon depinning up to a second, higher 
threshold force where a dynamical re-coupling transition 
occurs. Finally, these authors also observe a sharp in- 
crease in the depinning threshold at the crossover or tran- 
sition from continuous to hysteretic depinning, not unlike 
that shown in Fig. [2] A strong crossover from elastic 
to plastic with increasing disorder strength, with an as- 
sociated sharp rise of the depinning threshold, has also 
been seen in a variety of two dimensional simulations, 
such as those by Faleski, et alm^ Macroscopic hysteresis 
has not, however, been observed in these two-dimensional 
models. Our work suggests that mean field models with 
strong disorder tend to overestimate hysteresis. In mean 
field there is no range of correlation lengths and hystere- 
sis will always occur when the system is driven from a 
strongly pinned incoherent phase, where all degrees of 
freedom depin independently at the single particle de- 
pinning threshold. Upon depinning, the system acquires 
long range order and becomes therefore much stiffer, so 
that when the force is ramped down it can remain in the 
sliding state down to much lower values of the driving 
force. 

Early transport experiments on current-driven vortices 
in NbSca showed S-shaped IV characteristics at high 
magnetic fields with a peak in the differential resistance 
as a function of driving current^ Other puzzling effects 
were observed in the region of the peak, including unusual 
frequency dependence of the ac response and fingerprint 
phenomena. These experimental findings were originally 
interpreted in terms of plastic depinning of the vortex 
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system and macroscopic coexistence of disordered and 
ordered bulk vortex phases. This interpretation was cor- 
roborated by a number of simulations in two-dimensions, 
where the crossover from elastic to plastic depinning is 
clearly seen as a function of disorder strength. For strong 
disorder the system exists in a disorder static phase that 
depins plastically and then undergoes a dynamical or- 
dering transition to a moving ordered phase. The peak 
in the differential resistance corresponds to such a dy- 
namical ordering transition and in simulations is clearly 
associated with a sharp drop in the number of topologi- 
cal defects in the driven lattice. More recent experiments 
have suggested, however, that the disordered phase is a 
metastable phase that is injected at the sample's edges 
and then anneals into the stable elastic phase as it gets 
driven into the sample i^Si^SiS This interpretation has 
been confirmed by comparing transport experiments in 
the conventional strip geometry, where the edge effect is 
always present, to experiments in a Corbino disk geom- 
etry, where the vortices are driven to move in concen- 
tric circular orbits in a disk-shaped sample, eliminating 
boundary effects. Although there is mounting experi- 
mental evidence that these edge contamination effects 
may indeed control much of the vortex dynamics ob- 
served in experiments, the comparison with simulations, 
where coexistence of bulk ordered and disordered phases 
is routinely observed, remains puzzling. Of course one 
important difference is that most of the simulations are 
carried out at zero or very low temperature, where the 
disordered phase may be artificially stabilized. 

Substantial phase slip effects have also been observed 
in CDW systems, especially at the contacts'*^, and have 
been associated with the "switching" observed in certain 
materials. The reported correlation between broadband 
noise and macroscopic velocity inhomogeneities also sup- 
ports the idea that in these systems the dynamics may be 
dominated by large scale plasticityi^ While the switching 
itself has also been explained as arising from the presence 
of normal carriers^®, phase slips seem crucial to account 
for the correlation between broadband noise and macro- 
scopic velocity inhomogeneities. 

Finally, similar behavior has also been observed in col- 
loids driven over a disordered substrate. Pertsinidis and 
Ling'"'^ have studied experimentally single layers of two- 
dimensional colloid crystal driven by an electric field over 
a disordered substrate. They observe plastic-like or fila- 
mentary flow of the colloids, with a velocity- force curve 
that is always convex upward and shows no hysteresis. 
Langevin simulations by Reichhardt and Olson^^ find a 
sharp crossover from elastic to plastic depinning as the 
strength of substrate is increased. Though the direct ap- 
plicability of our mean field model and results to experi- 
mental systems remains to be proven, this work lays out 
a detailed foundation for understanding the role of phase 
slips and topological defects on the dynamics of driven 
disordered systems. Preliminary numerical studies of the 
phase slip model in three dimensions, with a sinusoidal 
pinning potential, suggest that the depinning transition 



may not be hysteretic in the thermodynamic limit. This 
is similar to what suggested by studying the mean field 
with a broad distribution of pinning strengths, as shown 
in Fig. 1201 where the distribution of pinning forces the 
incoherent static (IS) phase. Clearly more work is needed 
to establish if such a finding is generic in finite dimen- 
sions. An important open question is whether the tran- 
sition from elastic to plastic depinning (with or with- 
out macroscopic hysteresis) is a crossover or is associated 
with some type of tricritical point, as suggested by the 
present and other mean field models. 

This work was supported in part by NSF grants DMR- 
9730678, DMR-0109164 and DMR-0305407. 

APPENDIX A: COHERENCE AT F = 

In this appendix we describe the calculation of the co- 
herence r(/i) of static states at F = 0. First we derive an 
expression for the function f{u) defined in Eq. Q for an 
arbitrary pinning force, Y{d). Once f{u) is known, the 
coherence is then obtained by solving the self-consistency 
condition, r = f{fJ.r). The calculation is complicated by 
multivalued solutions to the self-consistency equations, 
which leads to multiple metastable states. A consistent 
selection principle is applied, namely, choosing the coher- 
ence u to be maximal, given /x. The range of available 
metastable states is also used to determine fir, the value 
of coupling above which the depinning field is zero. 

1. Change of variables 

As discussed in Section Fill CI it is convenient to per- 
form a change of variables in Eq. ^ and integrate over S 
rather than over the random phase p. The function f{u) 
is then given by 

^^"^ ^hf"^^ {%) + ' ^^^^ 

where u ~ [ir. Since YiS) is 27r periodic, the integra- 
tion in Eq. (|Aip can be carried out over any 27r inter- 
val. Here we choose the interval [— tt, tt]. The change 
of variable allows us to evaluate j(u) analytically as the 
force balance equation, Eq. ((HJ, while transcendental in 
d{f3,u), is simply a linear equation in the phase P{d,u). 
We can therefore immediately write the solution (3{S,u) 
of Eq. substitute it in Eq. HA1|I . and evaluate the 
integral to obtain f{u). As we will see below, the only 
difficulty in carrying out this program is that the phase 
d{f3,u) is generally a multivalued function of /?. There- 
fore care must be taken in selecting the portion of the 
curve that must be included in the integral. The choice 
is dictated by the requirement that the imaginary part 
of the self consistency condition, which now reads 

r dS (^) sm[d + P{S, u)] , (A2) 
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be satisfied, and that the phase S{(3) span a full 2tt inter- 
val in f3. 

For static solutions and F — the balance equation 
© can be written as 

sin{S + ^ . (A3) 
u 

Since | sin(^ + /3)| < 1, the right hand side of Eq. (|X3|) 
must also be bounded in magnitude by one. This means 
that all solutions to Eq. ljA3|l must satisfy 

{u)<5<S 

max 

(u) , (A4) 

where 6max(u) is defined by 

\Y{dm..)\--u, (A5) 

or Smax{u) = \Y^^{u)\, with denoting the inverse 
function. Note that if Y{5) is non-monotonic in the inter- 
val [— TT, tt], as it is for instance the case for the soft cubic 
potential shown in Fig. ^b), then for u > |l^(7r)| there 
are two possible values of in the range [0,7r]. 

In this case 6max is defined as the smallest of these two 
values. At the end of this Appendix we will discuss the 
relevance of the second solution and demonstrate that it 
corresponds an unstable state. 
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FIG. 22: The figure shows the behavior of the phase S{f3) 
for three values of u. The n = ±1 half-sections are dashed, 
while the n — section is solid. Curve (a) corresponds to 
u < a and is single valued. Curves (b) and (c) are both 
multivalued and correspond to (b) a < u < Y{tt/2) and (c) 
u > Y{tt/2). The section /3o ends at the points ±Smax, where 
the half sections (3 = =pl begin. For curve (b) these points lie 
within the portion of the curve that must be included in the 
integral to determine f{u). For curve (c) they lie outside. S, 
denotes the non-zero value of the phase at (3 = —n. 



2. Metastable states 

For every fixed value of u, there is in general an infinite 
set of solutions for the phase S in the range [—Smax, Smax] ■ 
The corresponding solutions for the phase /3 as a function 
of S can be enumerated by indexing them with an integer, 
n. They are given by 

I3n{5, u) = -5 + mr+ (-1)" sm-\Y{S)/u) , (A6) 

where we only consider values of the function sin^^(a;) in 
the range [— 7r/2, 7r/2]. Since the calculation of /3„((5, u) 
and d{f3,u) is carried out at fixed u, from here on we 
will simply omit the u dependence in the argument of 
these functions. The typical behavior of the phase (5 as a 
function of /3„, for — 2 < n < 2, is shown in Fig.[S| 

The integral in Eq. (|A1I) must span a full period (in (3) 
of the 6{f3) curve. As evident from Fig.[Sl this always cor- 
responds to a pair of consecutive even-odd sections. Here 
we choose to work with the n = section, and the upper 
and lower halves of the n = —1 and the n — 1 sections, 
respectively. This choice is equivalent, for instance, to 
that of the n = and the full n = —1 sections (or n = 
and n = I), but it has the advantage of being symmetric 
about the origin. The chosen portion of the S{P) curve is 
displayed in Fig. |^ for three different values of u. The 
figure shows how the phase becomes multivalued as u is 
increased. 

For u < a, with a the linear slope of the pinning force 
Y{6) at 5 = 0, the phase S is single- valued, as in curve 
(a) of Fig. 1221 In this case integrating over a full period 



in P is equivalent to integrating over the entire curve, 
consisting of the full n — central section (solid) and the 
two n = ±1 half sections (dashed). Making use of the 
symmetry of the integrand about 6 = 0, we obtain 

f{u) = (^^^cos[^ + /?_i(5,u)] 

+ dS cos[<5 + (3oiS, u)] . (A7) 

Upon substituting the expressions for l3-i{S) and l3o{S) 
from Eq. ljA6|l in Eq. I|A7|I . we obtain 

/(m) = - / dS y/l - [Y{6)/u]^ , u<a. (A8) 
Jo 

When 6{f3) is single valued, the integral in Eq. H1(J|) over 
the entire period gives zero, so that the imaginary part 
of the self consistency condition is satisfied. 

When u > a, the phase S is multivalued, as exempli- 
fied in cases (b) and (c) in Fig. |221 In this case one 
can no longer simply integrate over the full curve in 
the range /3 G [— tt, tt]. Rather, one must select a por- 
tion, of measure 27r in /?, that satisfies the imaginary 
part of the self consistency condition, Eq. ljA2|l . As dis- 
cussed in Sec. IIII CI we choose the portion of the curve 
corresponding to the metastable states that would be 
accessed by adiabatically increasing u from zero. For 
i/j = 0, this choice corresponds to the connected part of 
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the S{P) curve lying between f3 = — tt and tt. This choice 
is odd about the origin and therefore automatically satis- 
fies Eq. (jSU. The phase 6 now has two values at /3 = tt, 
5 = and (5 = <5*, which is defined implicitly as the 
nonzero root of the equation 



Y{S^) = usin((5*) 



(A9) 



The value 5, is the desired upper limit in the integration 
over 6 in Eq. (jAlj) . When a < u < F(7r/2), correspond- 
ing to curve (b) in Fig. [52 the root (5, is smaller than 
5max E^nd the portion of the curve to be included in the 
integrand spans the entire (3o{S) section (solid line) and 
those parts of the /3±i((5) half sections (dashed) that he 
within f3 = [—tt, tt]. For this range of u values we find 



f{u) = - 



ds v/1 - iY{s)/uy 



dS v/1 - iY{S)/u)^ , 

a<u< \Y{tt/2)\ . (AlO) 



At u = \Y{tt/2)\, (5, = Smax- For u > \Y{n/2)\, cor- 
responding to the situation illustrated in curve (c), S^ax 
exceeds (5, and the portion of the curve to be included in 
the integrand only spans that part of the f3o{S) section 
(solid) that lies in /3 G [— 7r,7r], as seen from Fig. [221 In 
this case we obtain 



fin) 



- f'dS^l- {Y{5)/uY M 
u > \Y{tt/2)\ . 



(All) 



The three equations, Eqs. (|X5|) . (|A10|I and (|A11() . give 
the function f{u) at all u for an arbitrary pinning force, 
Y{S). It can be shown that when Eq. (|A8|) is expanded 
for small u, the perturbative result, Eq. H19() . is recovered. 

For a piecewise linear pinning force, with Y{d) = ~aS 
for -TT < (5 < TT, the integrals in Eqs. HA8|1 . HA10|1 and 
(|Alip can be evaluated analytically, with the result 



u 



f{u) 



2a ' 

|^(^ + H+coS<5,), 



u < a, 

a < u < aTr/2, 
u < an/2, 

(A12) 

where (5* = (u/a) sin((5*). The coherence r is then deter- 
mined by the solution of r = fifJ-r). For u < a the equa- 
tion for the coherence is r = fir/fi^, where /x^ = l/(2a). 
If /z 7^ the only solution is r = 0. For /i = /i„ the 
equation is satisfied by any nonzero value of r consistent 
with u < a, or equivalently r < 1/2. Thus, at fi = fi^^ the 
coherence jumps discontinuously from zero to the value 
ro = 1/2. By expanding f{u) for u ^ a+ we find that 
for fj. > 

r - ro cx (/i - ^„)^ . (A13) 
The full solution r as a function of fi is shown in Fig. [S] 
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FIG. 23: The figure shows the phase 5 versus /3„ at « = 0.8, for 
n = 0,±1,±2, for the continuous pinning force of Fig. ^c), 
with a = 3V3/(27r) and c = -3\/3/(27r^). The upper and 
lower branches, lying outside the range [— 
stable, while the central branch is stable. 



We now return to the question of the existence of so- 
lutions (5(/9) outside the range [—Smax, Smax]- This is rel- 
evant for pinning forces Y(5) that are non-monotonic in 
the interval [— tt, tt]. For such pinning forces the Eq. (|A9|) 
has two nonvanishing solutions. The smallest of these 
two solutions, J* defines the range of phases that have 
been used in the calculation of the coherence described 
above. Denoting the largest of the two solutions by (5„, 
we note that for u > |i^(7r)| there will also be solutions 
for the phase S lying in the ranges [(5i,,7r] and [— tt, —5u\- 
Examples of such solutions are shown in Fig. [23 for the 
soft cubic pinning force. The solutions outside the range 
—5max l£ 6 < Smax SiTe the top and bottom branches in 
the figure. It can be shown that such solutions are always 
unstable, while the center branch is stable. This is easily 
seen by plotting the total force Ftot = —u sin{5+f3)+Y{5) 
acting on a domain versus the phase 5, for a fixed value (3. 
The stable solutions of the force balance equation are the 
zeros of Ftot{5) with a negative slope, so that they corre- 
spond to minima of the total potential. The zeros with a 
positive slope are maxima of the potential and therefore 
represent unstable solutions. Of the two zeros shown for 
instance in Fig. [^ for (3 = 7r/2, only the left solution, 
which lies in the range [—Smax, Smax] is stable, while the 
right one is outside this range and is unstable. Changing 
the value of (3 would simply shift the curve of Ftot along 
the S axis, with the stable root always remaining inside 

the interval [-Smax, Smax]- 
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FIG. 24: Plot of Ftot versus 5 for /? = 7r/2, u = 0.8 for the 
continuous pinning force of Fig.Q^c), with a = 3\/3/(27r) and 
c = — 3\/3/(27r^). The equation Ftot = has two solutions. 
The left solution, with negative slope is stable, while the right 
one, with positive slope, is unstable. 



3. Derivation of \mt 

The number of metastable static states available to the 
system plays an important role in determining the depin- 
ning threshold. In general the system can exist in a large 
number of static metastable states and the function 5(/3) 
becomes more multivalued as u increases, as shown in 
Fig. 1221 The number of metastable states is not, however, 
a monotonic function of u as only values of b lying in the 
interval [— tt, tt] are acceptable solutions. The number 
of available metastable states increases with increasing 
u until (5maa;(it) = TT, Corresponding to w = lF(7r)|. As 
u is increased beyond |y(7r)| the number of metastable 
states decreases. When an infinitesimal force is applied, 
all the phases are pushed forward and an infinitesimal 
number of static metastable states becomes unstable as 
they can no longer satisfy the self-consistency condition. 
The system remains, however, pinned, provided there ex- 
ist other static states that are still metastable. When 
b^(u) — TT, the situation changes as there is only one 
metastable static solution that becomes unstable as soon 
as an infinitesimal driving force is applied to the system. 
The system depins as soon as > 0, i.e., the threshold 
force for depinning is zero. 

It can be seen from Eq. (|A9|) defining (5* that for pin- 
ning forces with |i^(7r)| > 0, (5* < tt for any finite u. In 
this case b„ approaches tt only in the limit u — > oo. Since 
r is always finite, it is only in the limit of infinite [i that 
the system approaches a perfectly ordered floating state 
and the depinning threshold force goes to zero. For con- 
tinuous pinning forces with ^(Tr) = 0, (5* = tt at a finite 
value of M = ut = Y'{ti). For u > ut, the system has 



only a single, albeit partially disordered, state available. 
This state becomes unstable upon application of an in- 
finitesimal driving force, and the system begins to slide. 
In other words, the threshold for depinning vanishes for 
all u > Ut or, equivalently, all fi > fix — UT/fiux)- 
Using Eq. (|A11|I we find the value of fix displayed in 
Eq. 

APPENDIX B: DEPINNING FORCE F^(fi) 

In this appendix we calculate the depinning force 
for hard and soft cubic pinning forces, of the type 
sketched in Fig. ^ These forces are given by Eq. 
with c > for the hard cubic force and c < for the 
soft cubic force. Due to the periodicity of the problem, 
we can restrict ourselves to any interval of i5 of range 27r. 
For simplicity we choose again 6 to lie in the [— tt, tt] in- 
terval. In this interval the force balance equation, with 
i/j = is 

= F-u sm{6 + (3)+ Y{6) , (Bl) 

and only solutions to Eq. (|B1|1 which satisfy — tt < 
6{P,u,F) < TT should be considered. 

As for the case F = 0, the transcendental nature of the 
force balance equation, Eq. (jBl|l . can be circumvented 
by integrating over 6 rather than over the phase /3 in the 
self-consistency conditions. Solving for (3{S, u, F) gives 
an infinite set of of solutions, labeled by an integer n, 

MS) = -S + mi + (-1)" sin-^ ( ^^ V ^ ) ' ^^^^ 
where 5 is restricted to lie in the range 

5min{u, F) < 5 < 5max{u, F) , (B3) 

with 

(5™„(u,F) = -y-i(F-7i), 

5^a.{u.F) = -Y-\F + u) . (B4) 

The solution must satisfy the real and imaginary parts 
of the self consistency condition, given by, 

r = f{u,F), (B5) 

0-^^ rf'5(^|^^sin[<5 + /3(<5,u,F)] , (B6) 

with 

dS cos[S + m F)] . (B7) 

Throughout the analysis we will be considering 5{[3^ u, F) 
for fixed values of u and F. We will therefore write 6 — 
S{p), with the dependence on u and F implied. 

As in the case F = 0, the phase S is generally a mul- 
tivalued function of p (see Fig. I25|l . We consider only 
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n = 




FIG. 25: The phase as a function of /3„ for -2 < n < 2 for 
two sets of values of (u, F), corresponding to single- valued (for 
(w = 0.1, F = 0.2)) and multi-valued (for {u = 0.35, F = 0.2)) 
solutions. Even n branches are drawn as solid lines and odd 
n branches are dashed. 



the metastable state corresponding to a connected por- 
tion of the curve 5{P) in the range S & [5l, Sr], and it is 
this portion that is integrated over in the self-consistency 
conditions. We focus in this particular state because it 
is the one that controls depinning. The points 6l and 
Sr bounding this portion are functions of u and F and 
may in general differ from Smm and Smax- They are de- 
termined by the requirement that the imaginary part of 
the self-consistency condition, Eq. (|B6|I . be satisfied and 
by the condition that the portion of the function d{l3) 
bounded by these points span a full 2tt interval in /3, i.e., 



P{dL) + 27r = (3iSR) 



(B8) 



The details on how the limits of integration are deter- 
mined and the corresponding portion of the solution for 
P{6) is chosen in each case are given below. 

After evaluating the coherence, we can then proceed to 
compute the force F| (/i) where the static coherent state 
becomes unstable and the system begins to slide. As F is 
increased at fixed u, the whole S{/3) curve shifts upward 
and both Smax and Sl increase. The number of static 
metastable states in the range S G [— tt, tt] decreases, until 
at the critical force F^ only one static metastable state 
remains. This occurs where the largest value of S on the 
connected portion, denoted by reaches tt, i.e., 



du{u,Ff) = TT 



(B9) 



Upon further increasing F the system depins. Equa- 
tion (jB9|l defines the boundary of stability of the coher- 
ent static state, i.e., the depinning threshold, and can 



be solved to obtain F^{u). It will be shown below that, 
depending of the value of u, the connected portion satis- 
fying the self consistency condition may or may not in- 
clude 5max- For small values of u it will and 5u — Smax- 
At larger values of u, the connected piece does not in- 
clude Smax and Su — S^. Finally, the depinning thresh- 
old F^{ii) as a function of /x is obtained by eliminating 
u between the equation for the coherence at threshold, 
r — f(u, F^) and the expression for F^{u) obtained from 
Eq. p9)l . 



1. F^{ii) for monotonic Y{S) 

The monotonic class consists of all hard cubic pinning 
forces and those soft cubic pinning forces which have |c| > 
o/tt^/^. Since the function Y{S) is monotonic, its inverse, 
y^^(x), is single-valued in the entire range of interest, 
-1 < a; < 1. 

A full period of S{P) corresponds to a pair of consec- 
utive even-odd sections in n. In Fig. (|25|l we show plots 
(with even sections shown as solid lines and odd sections 
shown dashed) of S versus /3„((5) for two pair of values 
{u,F), chosen so that in one case the solution is single- 
valued and in the other it is multi-valued. In both cases 
the curves lack the symmetry of those for F = 0. In 
general, the value of u at which S{P) becomes multival- 
ued depends on F. At this value, denoted by Usv{F), 
each odd PniS) develops an inflexion point at S = Sg. In 
particular, for n = I, this requires 



f df3i{S, u,,„F) 
\ dS 

d^piis, us^.FY 



dS^ 



, (BlOa) 
. (BlOb) 



Using Eq. (|B2() for f3i{S,u,F) we obtain the following 
pair of equations, 

(u,.)2 = [Y'{S,f+[Y"{5,)]\ (Blla) 
F = -Y{S,) - Y"{5,) , (Bllb) 

which can be solved to determine Ua^{F). 

For u < Usv{F) the function S{(3) is single valued, as 
shown in Fig. liGf a) . Integrating over a 27r interval of 
f3 is equivalent to integrating over a full odd and even 
section. We choose Sl = Smax, which requires 5r = Smax 
and automatically satisfies the imaginary part of the self 
consistency condition. The function /(w, F) is then given 
by 



ds(^^]cos[S + Po{S)] 
dS{^]cos[S + P,{S)].iBl2) 
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tion is then given by 







dS i ^ ] sm[5 + P^i{6)] 



dS ] sm[d + Po{S)] 

dd [ ^) sm[S + (3,(5)] , (B14) 



with the additional requirement 

(3.ii6L) + 2n^Pi{6R) . (B15) 

Once the values of Sl{u, F) and (5_r(u, F) have been 
obtained by solving Eqs. ljB14p and (|B15|) . the func- 
tion f{u,F), is computed using Eq. (|B5|) . which 
now has the explicit form 



FIG. 26; The phase 5 as a function of (3 for various values of 
u and F < F^{u). Also shown in each frame are the values of 
Smin and 5max defined in Eq. I|B4^ and the boundary points 
Sr and Sl of the connected portion of the function S{P) that 
is used to evaluate the integrals determining the coherence 
in each case. The four curves corresponds to the four cases 
discussed in the text: (a) u < Usv{F), where is single- 
valued. In this case we can choose Sl = Smax, which requires 
Sr = Smax- As F is increased, S^ax grows until S^ax = 5l = 
n at F = Ff. (b) Usv{F) < u < Ui{F); (c) ui{F) < u < 
U2{F); and (d) u > U2{F). 



Equation ljB12|l can be simplified as 



u < u,,{F) . (B13) 



For u > Usv{F), the function S{P) is multivalued. In 
this case there are multiple possible metastable states 
S{l3). We can use any one of these to calculate r{u,F) 
as long as the chosen state satisfies the imaginary part 
of the self-consistency condition, and lies in the range 
[— TT, tt], but as explained above we choose to focus on the 
one corresponding to a connected portion of S{P). As u 
is increased at fixed F, Sl{u,F) increases and Sji{u,F) 
decreases. For hard pinning forces Sn reaches Smin before 
Sl reaches S^ax- It is then convenient to distinguish three 
regions. 



1. Usv{F) < u < ui{F), where ui{F) is the value of u 
where Sr — Smin- In this region the connected por- 
tion includes all of the /3o {5) piece and some of both 
the (3-i{S) and Pi{S) pieces as shown in Fig. ECT b). 
The imaginary part of the self consistency condi- 



f{u,F) 



dS { ^ ) cos[S + p^iiS)] 
dS{^]cos[S + Po{S)] 



+ 



dS ] cos[S + Pi{S)] . (B16) 



2. ui{F) < u < U2{F), where U2{F) is the value of 
u where Sl = Smax- In this region the connected 
portion includes only parts of the (3q{S) and /3„i((5) 
pieces 

In this region Sl > Smax, but Sr < Srmn, as shown 
in Fig. I26r cl. The imaginary part of the self con- 
sistency condition is then given by 



(^^)sin[,5 + /3_i(5)] 
dS (^^ym[S + PoiS)] , (B17) 



where 



P.i{5l) + 2tt = /3o{Sr:) ■ (B18) 



This pair of equations yields Sl and Sr^ which can 
then be used to calculate f{u, F) as 



f{u,F) 



+ 



dS 



dS 



fd(3- 



V dS 

dpo 
dS 



cos[S + p-i{S)] 
cos[S + Po{S)] . (B19) 



3. u > U2{F). In this region the simply connected 
portion of the S{f3) curve only contains part of the 
n = branch, and none of the n = ±1 branches as 
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shown in Fig. l26r d). The imaginary part of the self 
consistency condition reads 

- dS (^^) sin[5 + MS)] , (B20) 

with 

PoiSL) + 27T = po{SR) , (B21) 
and the function f{u,F) is given by 

f{u,F) = J\s (^^yos[S + /3o{6)] . (B22) 

As discussed earher, the depinning force is defined 
by Eq. ljB9|) . i.e. it is given by the value of F where 
6u = TT. For all values of u < U2, we can obtain a 
simple analytical expression for F^ since in this region 



6u{u,F^) = dmaxiu,F^) — TT. Substituting in Eq. (|B4 



we obtain 

S.ma.{u,F^') = -Y-\Ff 

which is easily solved to give 



Ff{u) 



1 



U < U2{Ff 



(B23) 



(B24) 



For u > U2, Sjriax IS outsidc the connected portion of 
the curve included in the integration and Su — 5^- So 
threshold is reached when 6l = tt. In this case it is 
convenient to directly solve for the depinning threshold 
by setting 5l = tt and F = F^ in the self-consistency 
condition, which is given by 

0= / dS(——-^ 
with 



sin[5 + /3o(<5,u,Ff)] , 
(B25) 



PoiSL = tt) + 27r = Po{SB.iu, Ff )) . (B26) 

Together these two equations yield F^{u). In Fig. llTTI) we 
plot F^{u) vs u for the hard pinning potential, Y{x) = 
-{x + x^)/{tt + tt^). 

The method for obtaining F^(fj,) for monotonic soft 
pinning forces is analogous to that for the hard pinning 
force, except for one difference. In the case of a mono- 
tonic soft pinning force, the value of Sl reaches S^ax 
before Sr reached (5mm (the reverse takes place for mono- 
tonic hard pinning forces). This means that U2 < ui so 
that region 1 is now defined by Usv < u < U2, region 2 
by U2 < u < ui and region 3 by it > wi. Of course the 
single valued region remains u < Usv It is not difficult 
to see that only the expressions for region 2 will differ. 
In this region the imaginary part of the self consistency 
condition becomes 







dS 



dS 



d(3o 1 . 



d6 



sin((5 + /3i(^)) , (B27) 




FIG. 27: The force Ff{u), as a function of u for the hard 
pmnmg force, Y{x) = -(a;-|-a;^)/(7r-|-7r^). The arrows indicate 
the values ul^ , ul , and U2 separating the four different regions 
discussed in the text. These values are defined by the relations 
Wsu(-Ff«„)) = u*„, ui{F^{ul)) = ul and ii2(-Ff(M|)) = U2- 
The plot becomes nonlinear beyond 112 where the threshold 
goes from being determined by Smax = n to being determined 
by 5l = TT. 



which along with 



MSl) + 271^ 13,(6 r) 



(B28) 



determines Sl and S^;. The expression for f(u, F) in 
region 2 is now 

fiu, F) = dS (^^) cos[<5 + PoiS)] 

+ j^^'' d5 cos[5 + (3i{5)] . (B29) 



2. For Non-Monotonic Y{5) 

The method for obtaining F^{fi) for non-monotonic 
Y{5) is analogous to that outlined for monotonic F^{u). 
Matters are complicated, however, by the existence of 
additional unstable solutions of the kind discussed for 
_F = in Appendix A. In principle there is no difference 
in obtaining F^{fi); one must simply be careful to en- 
sure that only stable solutions are being considered. The 
differences in the calculation are quite technical and we 
spare the reader the details. 



APPENDIX C: NUMERICS 

To explore the phase diagrams of the mean-field model, 
we numerically integrated the equations of motion to de- 
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termine v and r as a function of F and /x. As seen in 
the main text and earlier appendices, the macroscopic 
behavior can depend on the preparation of the initial 
state. For N degrees of freedom i — 1,2, . . . , N, the (3i 
for most studies were set uniformly, Pi = ^i. We stud- 
ied several different initial conditions. One of the most 
frequently used was to set all 9i = (3i dX F = 0, which 
prepares the system in the incoherent static (IS) state, 
whenever it is stable. In order to prepare the system in 
a static coherent state, all phases would be set equal to 
zero. Coherent moving or static states were also prepared 
by starting from a high field F with, say, random initial 
positions 9i. (Incoherent moving states were prepared 
in some portion of the phase diagram. When preparing 
incoherent sliding states, we used = N degrees of 
freedom, with M distinct values for (3; the values of 9i 
for each (3 value were equally spaced in time based on the 
periodic single-particle (r — 0) solution to the equations 
of motion for the given [3.) Given the initial conditions, 
we typically computed v{F) and r[F) at fixed /i. This 
was done by integrating the equations of motion Eq. (jS)) 
using the fourth-order Runge-Kutta scheme. The force 
was raised in small discrete steps: after some amount of 
time icq at fixed force, v and r are measured and then 
F is increased (decreased) some small amount 5F. With 
this algorithm, the time average of the ramp rate dF/dt 
is given by 6F/tcq. In some cases, we fixed F and ramped 
fi up and down in a similar fashion. 

While the ramp rate and system size does affect the 
depinning force, the force at which v goes from zero to 
non-zero, we find generically that for ramp rates smaller 
than 10~^ and sizes TV greater than 256, we obtain results 
for both the incoherent and coherent depinning line that 
are relatively independent of actual ramp rate or sys- 
tem size and agree with analytical calculations. There is 
agreement even though the coherent depinning curve is 
analytically obtained using the assumption that u is adi- 
abatically increased. For the simulations, on the other 
hand, F or fi is increased (decreased) slowly. Adiabati- 
cally ramping /Lt is not necessarily equivalent to adiabat- 
ically ramping u since the former does not insure that r 
changes slowly, but we do find the correct coherent de- 
pinning line by sitting at a fixed F and ramping up /j,. 

The analytical analysis in Sees. V and VI provides us 
with the depinning line as approached from the pinned 
phase, but it does not give us insight into the nature of 
the depinning transition. For example, there could be 
hysteresis in v{F) or r{F) for cyclical histories in the 
force, for sufficiently large system sizes and arbitrarily 
small ramp rates. Hysteresis in the order parameters 
implies that the depinning transition is discontinuous. 
If there is hysteresis in v{F), then the depinning line 
as approached from the moving phase must be different 
from the depinning line computed in Sees. VI and V. To 
numerically search for hysteresis, we prepare the system 
in a coherent moving state and lower the force until the 
system stops. If this repinning line is different from the 
analytical depinning result, hysteresis between the static 
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FIG. 28: (a) Area of the hysteresis loop in v{F) as F is 
cycled from large values to zero and then back up again near 
fie ~ 0.75 for the scalloped potential. Different system sizes 
and ramp rates dF/dt are shown. Plot (b) is just a blowup of 
(a) very near 



(IS or CS) and moving phases is present and there is a 
region where the two phases coexist. 

For every potential investigated, we find that there is 
a range of < /j, < /ie where there is a coexistence of 
the moving and stationary solutions. In general, there is 
hysteresis between coherent moving (CM) and incoher- 
ent static (IS) phases. For the piecewise linear pinning 
force, the hysteresis extends into the coherent pinned CS 
region. In other words, the coherence r jumps from one 
finite value to another at the depinning transition and 
there is hysteresis in both r and v {F^ ^ F^). The nu- 
merical evidence for this is shown in Fig. [2H1 which shows 
the area of the hysteresis loop, dF [v^{F) — v\F)] , 
where (F) and (F) are the histories of v{F) for ramp- 
ing the field down or up, respectively. The amount of 
hysteresis, as measured by this quantity, is independent 
of system size and dF/dt, which suggests that the simu- 
lations are close the adiabatic and infinite- volume limit. 
There is a jump down in the area of the hysteresis loop 
when ^ exceeds fiu, but the area is still non-zero for 
^ > Ai„. 

For the hard potential, with the history described 
above, r jumps to zero when the system becomes pinned. 
When the drive is increased back up again, the system 
depins at a different F^ when fi < /ig. However, we do 
not observe hysteresis between the CM and CS phases. 
In fact, the hysteresis when ramping F vanishes suddenly 
at ^ = fie- See Fig. [2^1 This is because the slope of 
the coherent depinning line starts to increase rapidly at 
{fie,Fe) and eventually becomes infinite before curling 
over to possible hysteresis. Above the point at which the 
slope becomes infinite, the analytic calculations suggest 
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FIG. 29: Area of the hysteresis loop in v{F) near fi^ for the 
hard case, where He is the intersection of and F^. Different 
system sizes and ramp rates are shown. 



that coherent depinning can be observed by increasing 
fi at fixed F. This was verified numerically. For the 
soft-potential cases tested we did not observe hysteresis 
between the CS and CM phases. Hysteresis is only ob- 
served between the IS and CM phases. 



APPENDIX D: STABILITY OF THE IM PHASE 

In this section we investigate the existence of a sta- 
ble incoherent sliding (IM) phase. We note that the ve- 
locity of a single degree of freedom is always a periodic 
function of time. To obtain a constant steady state ve- 
locity for a collection of incoherent degrees of freedom, 
we assume that at some initial time ti the i-th degree 
of freedom is at the minimum of its own potential well, 
which in turn is randomly shifted by f3, and perform an 
average over the random starting times ti. These are 
chosen to be random variables uniformly distributed on 
the interval [0, P], with P being the period, that is the 
time over which the phases advance by 27r. This pro- 
cedure guarantees that we sample uniformly all possible 
incoherent moving states such that the system reaches a 
steady state. 

Proceeding as in the study of the incoherent static 
state, we assume that the IM phase exists and study 
its stability. The self-consistency condition that must be 
satisfied by the mean field solution is given by 



r{t)e 



p p 



" dp 
2^' 



(Dl) 



where the phase 9{t — ti, (3) is the phase at a time t > ti 
obtained by solving the equation of motion. In eq. (|D1|I . 
both the coherence r and the mean phase ip are functions 



of time. As in the static stability analysis, we let 5{t — 
ti, p) = 9{t -ti,(3)- (3. We then perturb the IM state at 
the time t — 0, with a perturbation of the form 



Sp{-t„ 13) = 5o{~t,) ~ e sin {(3 + 5o(-ti)) 



(D2) 



with e << 1. After inserting this in Eq. I|D1|I . we evaluate 
the right hand side aXt — Qio 0{e), obtaining r{t = 0) = 
e/2 and V'(i = 0) = 0. We then use this to compute f(0) 
to linear order, with the result 



^(0) = 2 



f + J^§Y'{So{-t.)) 



(D3) 



If r{t) is monotonic in time, then its stability is entirely 
determined by the sign of r(0). We would then con- 
clude that there is a critical value fJ,c{F) of the coupling 
strength below which the IM phase is stable, with 



Mc(F) = -| J^^Y'i6o{~t, 



)) 



(D4) 



With a change of variable from ti to 6q (using the equa- 
tion of motion), one finds that /ic = for all F > Fgp for 
any continuous pinning force (since Y(Tr) = 0). For dis- 
continuous pinning forces, however, we can evaluate the 
integral in Eq. I|D4(I by splitting the integral in a contri- 
bution from the smoothly varying part of Y'{So) on the 
interval [— tt, tt) and a jump at So = tt. This gives 



Mc(P) = -p 



ln( 



F+\Yi7r)\ 
F-\Y{7:)\ 



)-2 



y(^) 



F+\Yi-n)\ 



(D5) 



where P is a function of F. For the piecewise linear force, 
one can evaluate P and find fJ,c{F) > for some F > Pgp- 
The critical value of ndF) is given by 



lJLc{F) = Ait, 



(D6) 



In the limit of large P, fic approaches zero. As F ap- 
proaches Pgp, on the other hand, /ic = fiu- In other 
words, the IM stability curve abruptly ends at (/iu,Psp) 
as there can be no IM phase for any F less than Pgp. 
A transition from an incoherent to a coherent moving 
phase was indeed obtained theoretically by Vinokur and 
Nattermann^- in a model of for layered charge density 
waves and also observed by Olson et al^^ in numerical 
simulations of layered superconductors. For strong dis- 
order, these authors found a transition as the drive is 
increased from a 2D state of decoupled moving layers to 
a 3D state where the moving layers become coupled. Our 
short time results suggest that a similar transition may 
occur in the isotropic system studied here. However, our 
numerical studies indicate that this transition may be an 
artifact of the short time analysis. When testing the sta- 
bility of a system prepared in the IM phase numerically, 
we find that r{t) is generally not a monotonic function of 
time. Furthermore, a perturbation of strength e always 
destabilizes the IM phase in the limit of large system 
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size, unless the strength of the perturbation is made to 
decrease with system size. Finally, we verified that the 
IM phase remains unstable if the somewhat artificial av- 
erage over the starting times U is replaced by an average 
over a narrow distribution of pinning strengths. Given 



these numerical findings, we conclude that the IM phase 
is typically unstable in the isotropic mean field model 
studied here, although of course we cannot rule out that 
the system could be prepared in such a state by some 
special initial condition. 
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